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PREFACE 


This paper is the nineteenth in a series of reports growing out of 
studies of radar cross-sections at the Willow Run Laboratories of the 
Engineering Research Institute of The University of Michigan. The 
primary aims of this program are: 


1. To show that radar cross-sections can be determined 
analytically. 


2. To elaborate means for computing cross-sections of 
various objects of military interest. 


3. To demonstrate that these theoretical cross-sections are 
in agreement with experimentally determined values. 


Intermediate objectives are: 


1, To compute the exact theoretical cross-sections of 
various simple bodies by solution of the appropriate 
boundary-value problems arising from electromagnetic 
theory. 


2. To examine the various approximations possible in this 
problem and to determine the limits of their validity and 
utility. 


3. To find means of combining the simple-body solutions in 
order to determine the cross-sections of composite bodies. 


4. To tabulate various formulas and functions necessary to 
enable such computations to be done quickly for arbitrary 
objects. 


5. To collect, summarize, and evaluate existing experimental 
data. | 


Titles of the papers already published or presently in process of publica- 
tion are listed on the preceding page. 


The major portion of the studies summarized in this paper were 
performed under Contract No. AF 04(645)-33, a prime contract of the 
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Western Development Division of the Air Research and Development 
Command (University of Michigan Project 2428); the contract was under 
the technical management of the Ramo-Wooldridge Corporation. Appen- 
dix F was done under Purchase Order No. $-96527 (University of 
Michigan Project 2355) sponsored by the Radio Division of the Bendix 
Aircraft Corporation. Other work was done under Purchase Order 

No. 7540 (University of Michigan Project 2360) under the sponsorship 

of the Ramo-Wooldridge Corporation. 


K. M. Siegel 


SECIRE TT 


S/ECIRIE 1 
UNIVERSITY OF MICHIGAN 
2428-3-T 


TABLE OF CONTENTS 


section : Title Page 
Studies in Radar Cross-Sections ii 
Preface lil 
List of Figures vi 
List of Tables Xi 
J Introduction and Summary ] 
Il The Radar Cross-Section of the Rudolph Missile 4 
Il Reduction of the Detectability of the Rudolph 
Missile 10 
IV Continuation of the Basic Analysis Involving the 
Cross-Sections of Ballistic Missiles 15 
V Conclusions and Recommendations for Future Work 16 
Appendix 
A The Radar Cross-Section of the R-W Re-entry 
Bodies (Configurations I — VI) 18 
B Effect of Perturbation to the Shape of the Re-entry 
Body on the Cross-Section in the 200 to 3000 Mc 
Range 58 
C Camouflage Materials 84 
Utilization of Fission Waste Products for Decoy 
Purposes 85 
1D Decoys 102 
F Further Analysis of The University of Michigan 
Drop-Test Data 106 
G The Cross-Section of the Tanks 161 
H Ionization and Related Topics 177 
References 184 
V 


Sle IRE 1 


Number 


2-1 
Le 


2-3 


2-4 
£25 


\oN) Ww 
' I 
roe 


ty — — 
i t t 
— ee) — 


2-2 
2-3 
2-4 
2-5 
2-6 
2-7 


pe ee ae ee ae oe 


SEC IRE 


UNIVERSITY OF MICHIGAN 


2428-3-T 


LIST OF FIGURES 


Title 





Typical Trajectories of an ICBM 


A Typical Parameterization for the Rudolph 
Warhead 


Cross-Section of One of the Parameterizations for 
the Rudolph Warhead 


The Terminal Propulsion Unit (Tanks) 


Cross-Section of the Tanks (Warhead Removed) 
at 200 and 270 Mc 


Configuration I with Perturbations 


Theoretical Nose-On Cross-Sections as a 
Function of Frequency 


Configurations I, II, and III 
Configurations IV, V, and VI 

Radar Cross-Section of Configuration I 
Radar Cross-Section of Configuration II 
Radar Cross-Section of Configuration III 
Radar Cross-Section of Configuration IV 
Radar Cross-Section of Configuration V 
Radar Cross-Section of Configuration VI 


Comparison of the Cross-Section of Configura- 
tion VI with the Cross-Section of a Sphere of 
Diameter Equal to the Maximum Diameter of 
Configuration VI 


R. R. D. E. Experimental Data 


SE CIR 1 


13 
19 
20 
26 
27 
28 
29 
30 
31 


32 
34 


SECRET 
UNIVERSITY OF MICHIGAN 
2428-3-T 


LIST OF FIGURES (Cont. ) 





Number Title Page 
A.2-9 Experimental Data — Sperry Gyroscope Company oe 
A. 3-1 Cross-Sections of Configurations I, II, Il], and A 

at 3000 Mc (Method I) 36 
A. 3-2 Cross-Sections of Configurations I, II, II, and A 

at 3000 Mc (Method II) 37 
A. 3-3 Cross-Sections of Configurations I, II, II], and A 

at 1000 Mc (Method I) 38 
A. 3-4 Cross-Sections of Configurations I, II, IJ, and A 

at 1000 Mc (Method II) 39 
A. 3-5 Cross-Sections of Configurations IV, V, VI, and 

B at 3000 Mc (Method I) 40 
A. 3-6 Cross-Sections of Configurations IV, V, VI, and 

B at 3000 Mc (Method II) 4] 
A. 3-7 Cross-Sections of Configurations IV, V, VI, and 

B at 1000 Mc (Method J) 42 
A.3-8 Cross-Sections of Configurations IV, V, VI, and 

B at 1000 Mc (Method II) 43 
Ad =9 Cross-Sections of Configurations I — VI 44 
A.4-]1 Carrot and Cone Models 46 
A.4-2 Cross-Section vs. Aspect Angle for Carrot and 

Cone Data 47 
A.4-3 Largest, Smallest, and Average Cross- 

sections — Carrot and Cone Experimental Data 48 
A.4-4 7-OC Type Models 50 
A.4-5 Largest, Smallest, and Average Cross- 

Sections — 7-OC Experimental Data 52 

vii 


5[— IRE 


S/[ECIRE 1 


UNIVERSITY OF MICHIGAN 
2428-3-T 


LIST OF FIGURES (Cont. ) 





Number Title Page 
A. 4-6 Cross-Section Data at 200 Mc (Vertical 

Polarization) 54 
A.4-7 Cross-Section Data at 200 Mc (Horizontal 

Polarization) — I 55 
A.4-8 Cross-Section Data at 200 Mc (Horizontal 

Polarization) — II 56 
B. 1-1 Configurations A and B 59 
B. 1-2 Configurations C and D 59 
B.1-3 Configurations E and F 60 
B. 1-4 Configurations G and H 60 
B.1-5 The Rounded Step in Configuration G 61 
B. 2-1 The R-W Re-entry Designs and Configurations 

A and B 63 
B.2-2 Cross-Sections of Configurations A — H at 

3000 Mc 72 
B. 2-3 Cross-Sections of Configurations A — H at 

1000 Mc 73 
B. 3-1 Nose-On Cross-Sections of Configurations I — VI 

and Configurations A — F at 200 Mc 76 
B.3-2 Nose-On Cross-Sections of Configurations I — VI 

and Configurations A— F at 350 Mc 77 
B. 3-3 Nose-On Cross-Sections of Configurations I — VI 

and Configurations A — F at 500 Mc 78 
B. 3-4 Cross-Section of Ogives at 200 Mc 80 

Vili 


SE CIE 1 


Number 


4-1 
4-2 
4-3 
4-4 
4-5 
4-6 
rey 
4-8 
4-9 
4-10 
4-1] 
4-12 
4-13 
4-14 
4-15 
4-16 


.4-17 


SECRET 
UNIVERSITY OF MICHIGAN 
2428-3-T 


LIST OF FIGURES (Cont. ) 


Title 





Cross-Section of Configuration I Compared with 
the Cross-Sections of Configurations G and H at 
200 Mc 


Theoretical Nose-On Cross-Sections as a Function 
of Frequency for Configurations A, C, E, andG 
Compared with the Cross-Sections of Configura- 
tion I 


S() for Object 59 S-Band 
S(w) for Object 59 S-Band 
S(w) for Object 59 X-Band 
S(¥) for Object 59 X-Band 
S(w) for Object 61,S-Band 
S(w) for Object 61 S-Band 
S(v) for Object 61 X-Band 
S(w) for Object 61 X-Band 
S(w) for Object 62 S-Band 
S(w) for Object 62 S-Band 
S(w) for Object 62 X-Band 
S(w) for Object 62 X-Band 
S(w) for Object 64 S-Band 
S(w) for Object 64 S-Band 
S(w) for Object 64 X-Band 
S(w) for Object 64 X-Band 
S(w) for Object 65 S-Band 


SE (IRE 1 


Page 


82 


83 
137 
137 
138 
138 
139 
139 
140 
140 
141 
141 
142 
142 
143 
143 
144 
144 
145 


SEC IRE 1 


UNIVERSITY OF MICHIGAN 
2428-3-T 


LIST OF FIGURES (Cont. ) 





Number Title Page 
F.4-18 S(w) for Object 65 S-Band 145 
F.4-19 S(w) for Object 65 X-Band 146 
F. 4-20 S(w) for Object 65 X-Band 146 
G. 1-1 Geometry of ''Tanks" Used in Computations 162 
G.2-1 A 5/6-Scale Model of the 7-OC Booster 163 
G.2-2 Theoretical Cross-Sections of the Tanks 165 
G.2-3 Cross-Section of the Tanks at 90 Mc 166 
G.2-4 Cross-Section of the Tanks at 270 Mc 167 
G. 3-1 Cross-Section of the Tanks at 200 Mc (Method 2) 169 
G.4-1 Geometry of the Modified Tanks Used for Com- 
putational Purposes 174 
G.4-2 Cross-Section of the Tanks at 200 Mc (Method 3) 176 


SEGRE I 


Number 


A, 


2-1 


A.4-1 


Qn ate ye ye yee ye tse et yt wD 


3-1 
3-2 
2-1 
eee 
.2-3 
.2-4 
2-5 


2-6 


.a-7 
2-8 
2-9 
.2-10 


4-] 
4-2 
4-3 
4-4 


S|ECIRE 
UNIVERSITY OF MICHIGAN 
2428-3-T 


LIST OF TABLES 


Title 





Nose-On Cross-Sections of Configurations I — VI 


Cross-Section Data for Configurations Shown in 
Figure A.4-4 


Nose-On Cross-Sections of Configurations A —F 
Comparison of Dimensions 

Frequency Distribution — Object 59 (S-Band) 
Frequency Distribution — Object 59 (X-Band) 
Frequency Distribution — Object 61 (S-Band) 
Frequency Distribution — Object 61 (X-Band) 
Frequency Distribution — Object 62 (S-Band) 
Frequency Distribution — Object 62 (X-Band) 
Frequency Distribution — Object 64 (S-Band) 
Frequency Distribution — Object 64 (X-Band) 
Frequency Distribution — Object 65 (S-Band) 
Frequency Distribution — Object 65 (X-Band) 
R(nT,) — S-Band — Swept Back Fins 

R(nT o>) — X-Band — Swept Back Fins 

R(nT,) — S-Band — Square Fins 

R(nTo) — X-Band — Square Fins 


Comparison of Dimensions Between the Tanks 
and the 5/6-Scale Model of the 7-OC Booster 


Page 


24 


116 
117 
118 
127 
129 
131 
134 


163 


SECRIEW 


UNIVERSITY OF MICHIGAN 
2428-3-T 


INTRODUCTION AND SUMMARY 


This report is a continuation of the radar cross-section study of 
ballistic missiles begun in Reference 1. In this report the radar cross- 
sections of a new ICBM for several parameterizations have been found 
for aspect angles out to 60 degrees off-nose, from 200 Mc through 5- 
band. It is concluded that there are several methods, which appear to 
be feasible from the cross-section point of view, for significantly re- 
ducing the radar cross-section of this new warhead; future systems 
analysis will determine if these changes will unduly disturb the over-all 
system. A major point established in this report is that fission waste- 
products are not suitable for camouflage purposes. 


The new ICBM warhead referred to above was designed as the re- 
sult of studies carried out by leading experts in hydrodynamics, ther- 
modynamics, aerodynamics, and upper-atmosphere physics under the 
leadership of the Ramo-Wooldridge Corporation. They felt that the 
third stage of anICBM as previously designed might vaporize upon 
re-entry into the atmosphere. Asa result, new parameterizations of 
the third stage of the ICBM were developed by the staff members and 
consultants of the Ramo-Wooldridge Corporation. The University was 
requested to determine the radar cross-section of these vehicles. The 
cross-sections determined and the analyses performed are presented 
in Appendix A; a brief summary of the results is presented in the body 
of this report (Sec. II). 


Enough analysis of the countermeasure problem has been performed 
to indicate that, if the modifications suggested in this report are made, 
the radar cross-section of the new ICBM warhead will be significantly 
reduced. The modifications involve certain changes in the shape of the 
vehicle while it is out of the atmosphere, without disturbing the vehicle's 
geometry upon re-entry into the sensible atmosphere. It is assumed that, 
for effective warning,the search radars must look for the ICBM long be- 
fore re-entry takes place. 
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Design of perturbations to Rudolph! to perform this function has been 
carried out; scale drawings of the new shapes and perturbations to them 
and computations of their cross-sections appear in Appendix B. It is shown 
there that perturbing the Rudolphshape inthe proper manner canresult in 
a reduction of the nose-on cross-section by a factor of 10° at S-band. 


The perturbation method of reducing the detectability of Rudolph 
requires that the warhead be in approximately a nose-on attitude to the 
enemy search radar. This in turn requires knowledge of the location 
of the enemy search radars and means to point the warhead in the proper 
direction while it is in flight out of the atmosphere. There are three 
other methods of reducing the detectability of an ICBM which are dis- 
cussed here: decoys, camouflage, and chaff, 


Since the aircraft decoy problem has been studied quite thoroughly 
(Ref. 2), little additional work seemed necessary from the reflectivity 
point of view to design decoys for the ICBM with the exception of con- 
sidering certain problems in the decoy field which are essentially dif- 
ferent for aircraft in the atmosphere and for missiles out of the atmos- 
phere. As aresult, only two efforts other than the perturbation- 
protuberance method were made in this field. These involved the analy- 
sis of the latest proposed camouflage material for aircraft and missiles 
(App. C), and the possibility of the use of radioactive wastes to serve 
either as an ionized medium itself or as an inducer of ionization in the 
atmosphere (App. D). Reference 2 contains a discussion of the possibil- 
ity of the defense determining the difference between decoys and aircraft. 
Since the same methods of differentiating decoys from aircraft are avail- 
able for use in defense against missiles, and since it was felt that these 
problems warrant further study, these methods are summarized in 
Appendix E. 


ISince no name had been given to the new third stage of the ICBM, 
the authors have coined the name Rudolph because of the "red nose’ 
upon re-entry. 
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One subject discussed in Appendix E in regard to decoys is the neces- 
sity for foreseeing the possible methods available for distinguishing 
between decoys and actual missiles or airplanes. One proposed method 
uses the differences in the fine structure of the radar return (scintilla- 
tion and glint) to make this distinction. Since the auto-correlation func- 
tions and power spectra associated with the radar cross-sections of bal- 
listic missiles have not been obtained previously, and since this infor- 
mation might prove important in more sophisticated methods of detection, 
these quantities were computed using The University of Michigan drop- 
test data. The drop-test experiment is summarized in Reference 3; the 
experimental cross-section data is presented in Reference 1; and the 
analysis of this data which yields the power spectra and the auto- 
correlation functions is presented in Appendix F. 


One must consider the cross-section and trajectory of the terminal 
propulsion unit (tanks) as compared with the cross-section and trajectory 
of the third stage to see if the knowledge of one will give away the loca- 
tion of the other. As a first step in this direction the radar cross-section 
of the tanks was determined. Also, some analysis was made of differ- 
ent ways of changing the enclosure which previously housed the third 
stage, in order to reduce the nose-on cross-section of the tanks. Com- 
putations for the tanks are given in Appendix G. 


Additional study has been carried out on the analysis of reflections 
from ionic trails such as those from meteors and, since it is clear that 
these trails are worthy of strong consideration by the defensive people 
in the ICBM field, a summary of thinking on this problem and on other 
related topics is given in Section IV and Appendix H. 


Recommendations for future research on the radar-reflection prob- 
lem for offensive ICBM study appear in Section V. 
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IT 


THE RADAR CROSS-SECTION OF THE RUDOLPH MISSILE 


The primary purpose in the design of the Rudolph warhead was to 
make sure that it would not vaporize upon re-entry into the atmosphere. 
It is expected that any low-drag supersonic missile traveling at the ve- 
locity of Rudolph would vaporize away uponre-entry. (Typicaltrajectories 
of an ICBM are shown in Figure 2-1.) The change in configuration was 
designed to produce a high-drag shape to replace the low-drag shapes 
previously under investigation (e.g., the 7-OC warhead). As a result, 
the small conical nose was replaced by a spherical one!, Of course, this 
aerodynamic difference changed the radar cross-section picture and at 
small wavelengths increased the cross-sections of the ICBM warhead. 
The detailed analysis of the cross-section work for the different para- 
meterizations of the Warhead is given in Appendix A. For the sake of 
clarity and brevity, one parameterization of the six considered in Appen- 
dix A is dealt with here; this particular third stage is referred to in App- 
endix A as Configuration I and is presented in Figure 2-2. 


In Appendix A the cross-sections of the various parameterizations of 
the Rudolph warhead are given out to 60 degrees off-nose for wavelengths 
between 10 cm and 150 cm (i.e., for frequencies between 3000 Mc and 
200 Mc). For these conditions Configuration I yields the smallest cross- 
sections. The cross-sections of Configuration I are shown in Figure 2-3. 


The cross-section of the terminal propulsion unit (the "tanks," Fig. 
2-4) of Rudolph is discussed in Appendix G. Several different variations 
of the surface of the enclosure exposed, when the Rudolph warhead has 


'The 7-OC warhead has a conical nose of half-angle 15 degrees. 
An earlier configuration, referred to in Reference 1 as the 'needle-nose 
warhead," had a conical nose with total-angle equal to 7 degrees. 
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41 min. flight time 


<@—— ]0° of latitude 





Surface of Earth 


Center of WA 
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FIG. 2-1 TYPICAL TRAJECTORIES OF AN ICBM _ (Information for this 
Figure was obtained from Ref. 4) 
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Rn = 0.15m 
Ry = 0.20m 
Ry =0.47m 
Ro =0.50m 





FIG. 2-2 A TYPICAL PARAMETERIZATION FOR THE RUDOLPH WARHEAD 
(CONFIGURATION 1) 


separated, are discussed. These variations are: 


a. The hole filled with an absorbing material so that only 
a negligible contribution to the cross-section is received 
from the interior of the hole. 


b. The hole lined on the inside with an inverted hemisphere. 


c. The hole considered empty up to the end of the hat section, 
so that for the nose-on aspect the return from the hole 
would be like that from a flat plate. 


d. The interior of the hole lined as the inverse of the rear 


portion of the warhead. 


The first and last of these variations yield the smallest nose-on cross- 
sections. 


The results (with the warhead removed) are summarized in Figure 
2-5, 
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FIG. 2-3. CROSS-SECTION OF ONE OF THE PARAMETERIZATIONS FOR 
THE RUDOLPH WARHEAD (CONFIGURATION 1!) AS A FUNCTION OF 
ASPECT AND FREQUENCY 
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REDUCTION OF THE DETECTABILITY OF THE RUDOLPH MISSILE 





There are four approaches to the problem of reducing the detectability 


of the Rudolph missile: 


l, 


If it is possible to choose the orientation of the warhead 
with respect to the radar within certain broad limits 

(5 to 20 degrees) when it is out of the atmosphere, the 
best method available seems to be to reduce the cross- 
section of the Rudolph warhead by the addition of 
structures which will vaporize readily upon re-entry 
into the atmosphere, thus re-establishing the desired 
shape. 


A large number of false Rudolphs could be produced so that 
the enemy could not discriminate between the actual warhead 
and the decoys. 


Camouflage materials could be applied to reduce the cross- 
section of the missile. 


The warhead could be hidden in an almost continuous column 
of ionized particles or in chaff. 


3.1 PERTURBATION OF THE SHAPE OF THE WARHEAD 


The discussion of perturbation of the shape of the warhead must be 
broken up into three cases: 


a. 


The assumption that the tanks would not have to be separated 
from the warhead until after the search phase was ended, 


The assumption that the tanks could be separated from the 


warhead by a large distance (the tanks could then be blown 
up), and 


The assumption that the separation would occur, but that 
the tanks would remain near the warhead. 
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Let us consider (c) first. In that case the defense could fix its search 
radars on the tanks and predict the trajectory of the warhead to within 
small bounds. Thus perturbation of the shape of the warhead by itself 
is not sufficient. 


Under case (b) the most economical method available for reducing 
the cross-section of the warhead is to put a thin sheet of tin around 
Rudolph (Fig. 3-1). It is expected that the tin will vaporize almost 
immediately upon re-entry, and as a result one will obtain the aerody- 
namic shape required for re-entry. If properly done, this will reduce 
the cross-section of the warhead while it is outside the atmosphere by 
a factor of approximately 10° at S-band in the nose-on region. In Fig- 
ure 3-2 we have plotted the average cross-section of the Rudolph war- 
head (Configuration I) within t5 degrees of nose-on as a function of 
wavelength, with and without the tin structural perturbation. Figure 
3-2 also contains the cross-sections of the 7-OC and Needle-Nose war- 
heads as a function of wavelength. In this figure observe that the cross- 
section of the Rudolph warhead has not only been decreased by a factor 
of approximately 10° at S-band, but, simultaneously, this cross-section 
has been reduced by a factor of about 10% at 200 Mc (App. B). 


Under case (a) the cross-section of the entire missile at the nose- 
on aspect will be comparable to the cross-section of the warhead alone 
above 200 Mc!. | 


3.2 THE USE OF DECOYS 


If additional load capacity is available in the Rudolph missile, this 
space can be used for decoys having the same cross-section as the Rudolph 
warhead, and adequate aerodynamic characteristics. For a discussion of 
the use of decoys see References 2, 5, and 6. This work is summarized 





1At lower frequencies the nose-on cross-section of the warhead 
alone will be less than that of the entire missile because the warhead is 
in the Rayleigh region. 
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FIG. 3-1 CONFIGURATION | WITH PERTURBATIONS 


The shaded portions indicate the perturbations which would vaporize 
upon re-entry. This perturbation is the one which is referred to as Configur- 
ation G in Appendix B. 
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in Appendix E. It is shown in Appendix D that fission waste-products 
cannot be used to enhance the cross-section of the decoys. 


3,3 THE APPLICATION OF CAMOUFLAGE MATERIALS 


The University of Michigan's examination of the camouflage field has 
consisted of a continuous but limited study, under Air Force Contract 
AF 33(602)-1070, over the past 18 months. A portion of this study is 
summarized in Reference 7. Additional studies were made since the 
completion of Reference 7. One of these additional studies consisted of 
an examination of thick and heavy materials (approximately 8 inches 
thick and 1 1/2 lbs/ft*) for which large reductions in cross-sections can 
be obtained for operation at frequencies above 500 Mc (Ref. 8). 


Recently a new camouflage material has been developed by Deutsche 
Magnesit, AG, Munich, Germany. This new material seems to be the 
best presently available for aircraft and missile applications. Tests on 
this camouflage material, carried out at the Rome Air Development 
Center, are discussed in Appendix C. This material is effective in re- 
ducing the back-scattering cross-section over small frequency bands. 
Unfortunately, its effectiveness may be due not so much to absorbent 
qualities as to the fact that it scatters the electromagnetic energy into 
angles other than the back-scattering direction. 


The above studies are still in process and it is too early to state 
that there is a camouflage material presently available which is suitable 
for use on missiles. 


3.4 THE USE OF CHAFF OR IONIZED PARTICLES 


It is clearly not possible to carry enough chaff to cover the space 
required to camouflage the Rudolph warhead over a Significant portion of 
of its trajectory. The alternative possiblity of using fission waste- 
products to hide the warhead can be discounted on the basis of the results 
contained in Appendix D. 
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IV 


CONTINUATION OF THE BASIC ANALYSIS INVOLVING THE 
CROSS-SECTIONS OF BALLISTIC MISSILES 


4,1 POWER SPECTRA ANALYSIS 


One method of discriminating between ballistic missiles and ballistic- 
missile decoys might be based upon scintillation and glint. Recently, for 
the Bendix Aviation Corporation, The University of Michigan computed 
the auto-correlation function and power spectra of the cross-section 
measurements obtained for one-third scale V-2 models (Ref. 1 and 3). 
The results of this investigation are given in Appendix F. 


4.2 IONIZATION AND SHOCK WAVE CONSIDERATIONS 


In Appendix H we have summarized the work presented in Reference 
1, and information on the analysis performed by the Massachusetts In- 
stitute of Technology onthe effects of ionization of the atmosphere by the 
missile. Also included is a discussion of the shape of shock waves formed 
by a missile as it leaves the atmosphere, as compared to the shape of the 
shock wave upon re-entry. 


4.3 EXPERIMENTAL WORK ON HIGH SPEED PELLETS 


It has been discovered recently by Dr. T. C. Poulter (of the Poulter 
Laboratories) that one may obtain experimental velocities in the neighbor- 
hood of 15 km/sec. These experiments were performed at Stanford Re- 
search Institute. The results are summarized in Appendix H. 
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V 
CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE WORK 


The radar cross-section study of ballistic missiles begun in Refer- 
ence 1 has been continued in the present report. In this report the radar 
cross-section of several parameterizations of the Rudolph warhead have 
been found from 200 Mc through S-band for aspect angles out to 60 de- 
grees off-nose. It has been concluded here that there are several methods 
for significantly reducing the radar cross-section of the Rudolph warhead; 
future systems analysis will determine if these changes will unduly dis- 
turb the over-all system. A major point established in this report is 
that fission waste-products should not be used for camouflage purposes. 


It is recommended that future work on the radar cross-section prob- 
lem for the ICBM should include: 


1, A thorough analysis of the decoy problem for the ICBM 
(similar to the thorough analysis given to the decoy 
problem for aircraft (Ref. 2)), 


2. A continuing analysis of the reflection characteristics of 
ionized media (as associated with problems of propulsion, 
as well as with ionization of the atmosphere during the 
re-entry portion of the missile trajectory), 


3 A continuation of the perturbation analysis presented in 
this report (The University of Michigan expects to do this 
on an extension of Air Force Contract AF 04(645)-33. ), 


4, An experimental check on the cross-section computations 
performed on the perturbations to Rudolph warhead shapes 
(It is planned also to include these experimental checks in 
the extension of Air Force Contract AF 04(645)-33. ), 


5. The determination of the bistatic cross-section of the 
ICBM (due to the advantages to the defense that bistatic 
operation would yield), and 
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6. Further analysis of camouflage materials, including 
the considerations of the idea of rotating the polari- 
zation of the reflected energy. 
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APPENDIX A 


THE RADAR CROSS-SECTION OF THE R-W RE-ENTRY BODIES 
(CONFIGURATIONS I —VI) 


A.1 INTRODUCTION 


This appendix contains the analysis involved, and the results 
obtained, in the computation of the radar cross-sections of the six Ramo- 
Wooldridge Corporation re-entry designs. The six configurations are 
shown in Figures A.1-1 and A.1-2.(Throughout this report, they are 
referred to as Configurations I-VI. ) 


Cross-section computations have been performed for 200, 350, 500, 
1000, and 3000 Mc for aspect angles (w) out to 60 degrees off-nose (in 
the case of 3000 Mc, the aspect angle » ranged from 0 to 180 degrees). 
Section A.2 contains the 200 —500 Mc data and analysis, and Section 
A.3 contains a discussion of the 1000 and 3000 Mc considerations. 


Section A.4 is devoted to a consideration of experimental data at 
200 Mc for shapes whose dimensions are comparable to those of Con- 
figurations I—VI. This is done because of the inapplicability of physical 
optics to the study of the cross-sections of these configurations at 200 
Mc. 


In estimating the cross-sections of the Configurations I—VI at 1000 
Mc and at 3000 Mc, the methods of physical optics were used exclusively. 
At the lower frequencies, due to the inapplicability of physical optics, 
other methods of approximation were required.’ For the nose-on aspect 
an approach based upon the creeping wave method was employed; in the 


lAt these lower frequencies the ratio of the characteristic dimension 
of the body to the wavelength is approximately equal to one. When this 
ratio is much greater than one, the methods of physical optics can be 
applied, and when the ratio is much less than one, the Rayleigh approx- 
imation is valid. 
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aspect range 0< ¥ <60°, cross-section estimates were obtained based 
upon the results of the study of experimental data (Sec. A.4), physical 
optics results, and the data obtained in the study of the nose-on case. 


The over-all results of any such approximation method should be 
accurate within a factor of 10. It is of course hoped that the approxima- 
tions will be better in many cases, but this maximum error is tolerable! 
because the range of detection varies as the fourth root of the cross- 
section. In other words, if the radar cross-section is known within a 
factor of 10, the range performance of the radar system is known within 
a factor of 1.78. 


A.2 CROSS-SECTIONS AT 200 TO 500 Mc 
A.2.1 Cross-Sections for the Nose-On Aspect (¥ = 0) 
cr A ES ED 


Each of the Configurations I-VI has a maximum diameter, D, and 
length, L, such that at 200 Mc, D/A = 0.67, and 0.93< L/A< 1.43; at 
350 Mc, D/d = 1.17, and 1.63< L/X< 2.50; and at 500 Mc, D/d = 1.67, 
and 2.33 L/X<3.57. Reference to the theoretical sphere curve in 
Kerr (Ref. 9) and the 10:1 prolate spheroid curve (Ref. 10), indicates 
that these values of D/\ and L/) are in the resqnance region. These 
ratios are of such magnitude that neither physical optics nor Rayleigh 
scattering can be considered to be applicable (although at 500 Mc we are 
approaching values of these ratios for which physical optics could be 
applied). Thus some other method of approximation must be applied to 
estimate the nose-on cross-sections of Configurations I-VI for these 
frequencies. 


The method applied is similar to the one which was successfully 
used by Franz and Deppermann (Ref. 11) for the cylinder and is the same 
as that applied in Reference 1 for the determination of the nose-on cross- 
section of the 7-OC warhead (The 7-OC calculations yielded results which 
agreed with experimental data to within a factor of 2.) Briefly, the 
method involves considering the back scattering cross-section to be the 
Sum of two contributions. One is due to the radiation scattered from 
the first Fresnel zone on the body or, in the case of a pointed body, 
from the tip, and the other is due to the radiation which creeps around 


'The accuracy of the results obtainable depends only upon the 
amount of time and money one wishes to expend in obtaining the results. 
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the rear of the body and comes back to interfere coherently with the 
radiation scattered from the forward part of the body. An estimate of 
the contribution from the rear was obtained in Reference 1 through the 
Study of sphere cross-sections. To do this the cross-section of a 
Sphere at a peak was written 


2 
o = a +A | ; 
peak g.0. rear 


where Cpeak is the sphere cross-section at a peak in the resonance 
region, Ago. is proportional to the field amplitude due to the first 
Fresnel zone, and Rane is the contribution due to the radiation which 
creeps about the rear of the body. Taking (A o.) =a’, where 

a = the radius of the sphere, and reading values of Speak Off existing 
Sphere graphs such as the one found in Kerr (Ref. 9, p. 453), it was 
found that in the region of interest 1 is given by 0. 033r¢ (d/a)!/2 : 


Applying this procedure to Configurations I—VI and considering the 
contribution from the rear to be approximated by that from a sphere of 
radius equal to 0.5 m, we obtain for each of the six configurations 


2 
(Avear) = 0.128 m* at 200 Mc, 
é 2 
(Ayear) = 0.0317 m*% at 350 Mc, and 
2 2 
(A cay) = 09-0130 m% at 500 Mc. 


From the physical optics methods applied for the 1000 and 3000 
Mc cases, the contribution from the first Fresnel zone is 


0.0707 m* for ConfigurationsI and IV, 
0.283 m* for Configurations II and V, and 
0.636 m* for Configurations III and VI, 
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for all three frequencies’ 


From these computed values ''average'' cross-sections can now be 
obtained by merely taking the sum of the two contributions, an upper 
bound by assuming complete reinforcement, and a lower bound by assum- 
ing complete interference between the two contributions; i.e., if o) is 
the contribution from the front and o> is the contribution from the rear, 
we have (7, + 72) as the "average'' value, (Vo - Yo) as the lower 
bound, and (Vo, 5 yo>)¢ as the upper bound. The results appear in 
Table A.2-1. 


At 500 Mc the ratios of L/\ are approaching values which would 
permit the application of physical optics. It is of interest to compare 
the values for 500 Mc given in Table A.2-1 with those which are obtained 
through the application of physical optics. The following average values 
are yielded hy the physical optics method for \ = 0.6 m: 


2 


Z 
ov - 0.17m : or = 0.39 ae CTI = 0.74 m ‘ 


Ty = 0.64 m*, oy = 0.84 m’, and oy; = 1.19 m*, 


It should be noted that these physical optics cross-sections are larger 
than those given in Table A.2-1 by factors ranging from one up to eight. 


A,2,2 Cross-Sections for the Off-Nose Aspects (0< ¥ < 60°) 
a ee ee 


The "creeping-wave" approach used for the nose-on aspect cannot 
conveniently be applied for ¥>0 since the use of a sphere to determine 
the contribution from the "rear'' for »>0 is a poor geometric approxima- 
tion to the actual shape. However, estimates of the manner in Which 
the cross-section varies over the interval 0 < » < 60° can be obtained 
through the use of 


1. The results obtained for \ = 0.1 m and } = 0.3m (Sec. A.3), 
2. The physical optics estimate at 500 Mc, 


'This contribution is given by rR’ where Ry equals the radius of 
the front spherical section of the missile (see Fig. A.1-1 and A.1-2). 
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3. The geometric optics prolate spheroid (a = L/2, b = D/2) 
curve in the region 0< ¥ < 60° (Ref. 10), and 


4. The results of the examination of experimental data 
(Sec. A.4) which indicated that at 200 Mc o(y) should 
not vary much in the interval 0< / <60° for these 
configurations, 


combined with the computed nose-on values given in Table A.2-1. Cross- 
section curves are given in Figures A.2-1 through A.2-6. Each figure 
contains, for a given configuration, the cross-section curves at 3000 Mc 
(A = 0.1 m) and 1000 Mc (A = 0.3m) (Sec. A.3), the curves for the prolate 
spheroid (geometric optics), and points obtained through physical optics 
at 500 Mc. With these as a guide, estimates of the cross-section curves 
at 200 Mc, 350 Mc, and 500 Mc are made with the o(0°) point determined 
from Table A.2-1. 


A.2.3 The Wavelength Dependence in the Region From 200 to 500 Mc 


Some idea of the manner in which the cross-sections of Configura- 
tions I—VI vary as a function of frequency in the region from 200 to 500 
Mc is obtained by an examination of Figures A.2-1 through A.2-6. How- 
ever, a better idea may be obtained by analogy with the variations with 
frequency of the cross-section of the theoretical sphere (Ref. 9) and the 
theoretical 10:1 prolate spheroid curves (Ref. 10). These shapes dis- 
play an oscillation of o as a function of \ in this frequency interval. It 
is to be expected that the nose-on cross-section of Configurations I-VI 
will oscillate as a function of \ in a somewhat similar manner. Figure 
A.2-7 shows the cross-section of a sphere (of diameter equal to the 
diameter of Configurations I—VI) as a function of 1/(2), together with 
the computed ‘average’, lower bound, and upper bound o at 200, 350, and 
500 Mc, for Configuration VI. Examination of this figure yields a fairly 
good picture of how the nose-on cross-section of Configuration VI varies 
in the 200 to500 Mc frequency region, which can be considered to be 
typical of the variations for all six configurations. 


Experimental data obtained by the Radar Research and Development 
Establishment (Ref. 12) and by the Sperry Gyroscope Company (Ref. 13) 
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add additional information relative to the oscillation of o with changing 
frequency. The RRDE data is reproduced in Figure A.2-8 and the 
Sperry data in Figure A.2-9. The configurations in these experiments 
have L/D ratios of from 6:1 to approximately 10:1. An examination of 
the above data, the theoretical sphere data, and the theoretical 10:1 
prolate spheroid data, indicates that Configurations I — VI have a total 
fluctuation in cross-section (as a function of frequency) of less than a 
factor of ten in the 200 to 500 Mc range. 


A.3 CROSS-SECTIONS AT 1000 AND 3000 Mc 


The cross-sections of Configurations I—VI at 1000 and 3000 Mc were 
computed by two methods. Method I yields "average cross-sections," 
and Method II yields cross-sections for which the oscillations due to 
phase differences are estimated. Since the computations were performed 
in connection with the perturbations to the six re-entry designs, the 
computational procedures and equations involved for both sets of cal- 
culations appear in Appendix B. 


The computational procedures followed closely the methods and 
techniques outlined in Appendix A of Reference 5, and the results obtained 
are shown in Figures A.3-1 through A.3-9. Figures A.3-1 through A.3-8 
cover the aspect interval 0< V< 60° and include the data for two of the 
perturbations (Configurations A and B) discussed in Appendix B. Figure 
A.3-9 contains 3000 Mc data in the aspect range from 60 to 180 degrees. 
Examination of the data contained in these figures indicates that of the 
six configurations, Configuration I has the smallest cross-sections at 
all of the frequencies and aspects considered. 


A.4 EXPERIMENTAL DATA AT 200 Mc FOR SHAPES COMPARABLE 
TO CONFIGURATIONS I-VI 


A.4.1 Introduction 


An examination of experimental data for configurations having body 
dimension to wavelength ratios similar to those for Configurations I—VI 
Should give information relative to the magnitudes of the Rudolph cross- 
sections. 
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FIG. A.2-8 R.R.D.E. EXPERIMENTAL DATA - EQUIVALENT ECHOING AREA OF 
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A.4.2 Experimental Data for Carrots and Cones 





The Microwave Radiation Company, Inc. has performed cross- 
section measurements for the Cornell Aeronautical Laboratory, Inc. 
under Contract DA-30-115-ORD-543 between the Cornell Aeronautical 
Laboratory and the Department of Defense. These experiments were 
conducted at a frequency of 9375 Mc on models which were 2, 4, 8, and 
32 wavelengths long. Two of the models involved were a ''Carrot' and a 
finite cone (Fig. A.4-1). This section contains some of the 2X model 
data scaled to 200 Mc. This results in full scale configurations which 
are all 3 meters long with a maximum diameter = 0.66 m and 0.75 m.’ 


To examine the effect of configuration shape on the model cross- 
sections, these two models are considered as four; i.e., the Carrot 
with the cone tip in front, the Carrot with the spherical bulb in front, 
the Cone with the tip in front, and the Cone with its base in front. These 
four configurations have approximately the same length-to-diameter ratio, 
but otherwise differ in many respects. 


The experimental data (Ref. 23) scaled to 200 Mc are displayed 
graphically in Figure A.4-2. 


Average curves can be fitted to this 200 Mc data which serve 
to predict most of the experimental points to within a factor of three. 
These curves, together with maximum and minimum curves, are shown 
in Figure A.4-3. For horizontal polarization the average curve is within 
a factor of three of all the experimental data for aspects out to 20 degrees 
from the front and is always within a factor of eight of the remaining 
data in the 25 to 60 degree aspect range. For vertical polarization the 
average curve is within a factor of three of all the data, except for the 
aspect of 40 degrees, at which the average is within a factor of five of 
the extremes. 


'This may be compared with lengths between 2. 16 and 1.40 meters 
and a maximum diameter of 1 meter for the six configurations of interest 
here. 
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A.4.3 Experimental Data for the 7-OC 


Static experiments on 7-OC-type missiles were performed by the 
Microwave Radiation Company, Inc. (MRC), and by the Evans Signal 
Laboratory (ESL) at the instigation of Project Wizard. These data 
(scaled to full 7-OC dimensions) are given in Reference 1. 


Using the 7-OC 75 Mc data appearing in Reference 1 and scaling it 
to 200 Mc yields data for the three configurations shown in Figure 
A.4-4, Here we have three configurations for which (at 200 Mc) 


L/xk = 1.54, 1.46, and 1.27; and 
D/d = 0.346, 0.360, and 0.346.’ 


As in the Carrot and Cone experiments, the number of configura- 
tions is increased (over the aspect range of interest) by considering 
the model rotated through 180 degrees to be a different configuration. 
Combining the data from the two sources (MRC and ESL) yields ten 
experimental curves both for horizontal and vertical polarization: 


1. 7-OC Model~—nose in front (MRC data), 

2. 7-OC Model—nose in front (ESL data), 

. 7-OC Model-—tail in front (MRC data), 

. 7-OC Model-—tail in front (ESL data), 

7-OC (Hemispherical Stern)—nose in front (MRC data), 
7-OC (Hemispherical Stern)—nose in front (ESL data), 
7-OC (Hemispherical Stern)—tail in front (MRC data), 


Oo ND OT _ W 


. 7-OC (Hemispherical Stern)—tail in front (ESL data), 


‘Compare with the 0.93 < L/A< 1.43 and D/\ = 0.67 ratios for the 
six R-W re-entry designs. 
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9. 7-OC (Flat Stern)—nose in front (MRC data), and 
10. 7-OC (Flat Stern)—tail in front (MRC data). 


Figure A.4-5 shows the largest, smallest, and average values 
obtained for the cross-sections of the ten shapes listed above by scaling 
the 75 Mc—7-OC data to 200 Mc. (Table A.4-1 contains the data from 
which Fig. A.4-5 was constructed. ) 


A.4.4 Summary of Experimental Data for Carrots, Cones, and the 7-OC 


By grouping the 7-OC and Carrot and Cone data and considering 
each of the models both from the front and the back, fourteen experimen- 
tal cross-section vs. aSpect curves are obtained for ten different con- 
figurations. In some cases these ten configurations differ markedly 
(compare Fig. A.4-1 and A.4-4), but the length (L) and maximum diam- 
eter (D) of each satisfy the conditions: 


1.27< L/\< 2.00 and 0.34< D/X<0.50. 


All of the data for vertical polarization are grouped in Figure A. 4-6, 
which shows the range of values of o (largest and smallest) as a function 
of aspect, together with an average curve. All of the data lie within a 
factor of five of the average curve. 


The horizontal polarization data for the five ''cone-nosed" config- 
urations are grouped in Figure A.4-7 ina similar manner. In this 
case the average curve fits the data to within a factor of six, except 
for the aspect angle range of 10 to 40 degrees. 


The remaining horizontal polarization data are grouped in Figure 
A.4-8. The average curve lies within a factor of five of all the data. 


The experimental data examined in this section indicates that the 
cross-section of a configuration does not vary to any appreciable extent 
with changes in configuration shape if the L/X and D/X ratios are 
restricted by the above inequalities, (This is especially true for vertical 
polarization.) The configurations involved in these experiments were 
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TABLE A. 4-1 


CROSS-SECTION DATA FOR CONFIGURATIONS SHOWN IN FIGURE A.4-4 AT 200 Mc (in m¢) 


Configuration source Pol. Aspect 
0° 10° 20° 30° 40° 50° ~~ 60° 


7-OC model—nose in front 








7-OC model~nose in front 

7-OC model—tail in front 

7-OC model—tail in front 

7-OC (Hemi. Stern)—nose in front 
7-OC (Hemi. Stern)—nose in front 
7-OC (Hemi. Stern)—tail in front 
7-OC (Hemi. Stern)—tail in front 
7-OC (Flat Stern)—nose in front 


7-OC (Flat Stern) —tail in front 
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all pointed at one end only. If the body were pointed at both ends, other 
data indicate that this relative constancy of the cross-section in the 0- 


to 60-degree aspect interval would be expected not to be as pronounced 
(see Ogive data in App. B). 


SEC IRE T 


UNIVERSITY OF MICHIGAN 
2428-3-T 


APPENDIX B 


EFFECT OF PERTURBATION TO THE SHAPE OF THE RE-ENTRY 
BODY ON THE CROSS-SECTION IN THE 200 TO 3000 Mc RANGE 
ee ee 


B.1 INTRODUCTION 


A part of The University of Michigan Project No. 2428 was the deter- 
mination of perturbations to the six R-W re-entry designs which might 
reduce the radar cross-section at 200, 1000, and 3000 Mc for aspect 
angles out to 60 degrees off-nose. 


In this perturbation study eight modifications, Configurations A 
through H, were considered. Configurations A, C, E, and G are mod- 
ifications of Configurations I-III, and Configurations B, D, F, and H 
are modifications of Configurations IV-VI. These eight perturbations 
are shown in Figures B.1-1 through B.1-4. 


Configurations A and B are obtained from the R-W re-entry designs 
by capping the nose with a cone; Configurations C and D enclose the 
re-entry shape except the rear hemisphere in an ogive; Configurations 
E and F are obtained from Configurations I-VI by capping the nose with 
a cone and filling in and rounding off the step; and Configurations G and 
H are obtained from the. R-W designs by capping both the nose and the 
tail with cones as well as filling in and rounding off the step. 


It is intended that these "caps'' be made of a material which will 
vaporize upon re-entry so that the configurations will reduce to the R-W 
design. Figure B.1-5 displays the process of rounding off the step in 
the case of Configuration G. 


In what follows the computational procedures are outlined and the 
results obtained are compared with the results for Configuration I. (It 
will be recalled from Appendix A that of the six R-W re-entry designs, 
Configuration I yields the smallest cross-section.) A summary of the 
nose-on cross-section data is given in Section B.4. 
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FIG. B.1-1 CONFIGURATIONS A AND B 
(a=20° for Conf. A and a=40° for Conf. B) 





FIG. B.1-2 CONFIGURATIONS C AND D 
(Ogives faired into the hemispherical stern) 
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FIG. B.1-3 CONFIGURATIONS E AND F 
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FIG. B.1-4 CONFIGURATIONS G AND H 
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FIG. B.1-5 THE ROUNDED STEP IN CONFIGURATION G 


B.2 COMPUTATIONAL PROCEDURES (1000-3000 Mc) 





The techniques employed in determining the cross-sections at 1000 
and 3000 Mc were based upon physical optics.’ As mentioned in Appen- 
dix A, the cross-sections of Configurations I-VI were determined by 
two methods, one yielding average cross-sections and the other taking 
into account phase differences. The same approach was used for Con- 
figurations A and B. 


'The basic theory involved is discussed in detail in Appendix A of 
Reference 5. 
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B.2.1 Computational Procedures for Configurations I-VI, A and B 


The formulas used in the computation of the cross-sections of Con- 
figurations I-VI, A,and Bare listed below by method.! (Since the pro- 
cedure was applied for aspect angles out to ¥ = 180 degrees for 3000 Mc, 
under University of Michigan Project 2360, the required formulas for 
the additional values of y are also included.) The meanings of the var - 
ious parameters involved in these formulas are defined in Figure B.2-1. 


B.2.1.1 Method I (0<¢ y < 1/2) 
For y = 0° 


_ 2 
o =7wR +7R, tan Qn: 


For 0°< pcan 





2 2 2 2 
g = TR a ec tan (a. + wv) +tan (a - | 
TT 
For aS p ~o. a 
2 R2 Z 
Gg = R 
_ n 8rsiny ae ay +) 


T 
For wp a) - a 


: Str fR,?/? - R,/eos?/*(a.) é 





og = 
9X sin a, cosa,, 
T T 
For = - ak p< 5 
, - >Re tan? (ay + ¥) 


87 Sin 


‘An improvement on the approach used here to estimate the "finite 
cone" contributions will be discussed in the next report on this contract. 
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FIG. B.2-1 THE R-W RE-ENTRY DESIGNS (CONFIGURATIONS | - VI) 
AND CONFIGURATIONS A AND B 
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B.2.1.2 Method II (0< y<1/2) 
For w= 0° 


gc=o0,+o.+2,/0,¢ cos ($, - ¢ 


] 2 1 2 2 


For 0<Y¥¢a 
n 


Pee, to, to, +2 /o,c, cos (¢, ~%.) 


+2 V1" cos (¢, - o,)+2 V°%3%% cos (4, -¢,) 


T 
Fer ae 


=o to, +2 ou, aye 
CSin, ie, vo, cos (¢, @): 


where 


2 2 
7, =7R, tan a P 


_AR>2 tan” (an + p) 

3 81 sin yp : 
XR, tan” (a, - ») 
Sie 


o e 
4 8 Sin w 
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$= 2 E: - i+ 
sina 
47 R, /cosy 
Oe Ae = (— -ij+— , 
l Xr sina, 2 


4n R, T 
2 \Xtana 2 ’ 
n 














¢ 


4nr R>cos (a, + v) T 
pe Dean Cee +]: and 


3 sin a, 


_ Ar R, COS (a, a y) 


I A 


4 ‘ 
dX sin a4 


B.2.1.3 Computations For + <WeTt 


In applying Methods I or II in the aspect range out to » = 180 degrees, 
the formulas listed below are employed. 


For =< pe ase 
2 2 t 


= - 2 V 
o oe ee Oe o% 
l 
Xr y®,8, [tan (v ~ a,) * tan 


| fae 
cos eB (RR, sing | : 


4m siny tany 
where 


XR, tan” (p - a) 
4 87 siny 
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2 
wT. = AR /8r sind tan » , and 
2 
7) = \R,/87 sing tan p . 
For » = 1/2 + ay 
3/2 
Str IR, -R, 


2 
3/2 Sage a, 
C cae. 
9 \siné a, COS a, 


Tw 
For > + se - a, 


o = Oo Py Se re 


+2 f°7"%4 cos (9. =O.) 
+2 Yoo. cos (o_ - $,) 
+2 y°7% cos (¢ =) 
+2 70,5, cos (¢, - #,) 
+2 74% cos (¢, - %) 
+2 V7. %, cos (¢,. - $) : 


where 
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TT 
zta<v<t ma (cont. ) 


2 
AR, tan” (¥ - a4) 
4 Sisiny 


2 

c= AR, /80 siny tan y , 
2 

= AR, /8 siny tan y , 


a 


4n Ry cosy 
o_ = oe) ae 
7 X Sin a, 2 


Ar R] cos O+ = ] T 
¢,=—— eae a ’ 
4 NN Sin a, 4 








5 = e, +wT , and 

4n(R, -R,)siny 
Oy ip. ae ee 
6 4 r 


Pore -aqey<t 


= ie t Pe oa t tt 
o vo dt ed ye aie 


+2 V y - ¢. 
v.07, COs (¢, ¢,) +2 v0, cos (¢., e.) 

rage) t! _ all t _ al 
V°775 cos (9. >.) +2 70, COs (9. by} 


+2 V ; - 9! V - 
To, cos (¢_ oy) +2 To. cos ($_ o.) 
+2 ; _ at " al 
V4" cos (¢, ?.) +2 V74"5 cos(¢ , o.) 
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T - a<ydm (cont. ) 


where 


+2 To, cos (64 - 6) +2 oye cos (¢, - ¢/) 
+2 "4% cos(¢, - ?,) +2 Vrs cos (o, = o:) 
ree, cos (¢, - 4) +2 Varo, cos (o, - ¢/) 
+2 Vow, cos (9, - #_) + 2V/oneg cos (¢,. ~ /) 
+2 |/o7o' cos (¢" “#6)+2 ore, cos ($7 - $4) 
+2 Vays, cos (44 “W)+ 2), cos (#) - $4) 
+2 Vopr, cos ley $9) ; 


2 
=w7R 5 


q t 


AR, tan@ (y - a) 


4 81 sin p ‘ 


2 
Feg!! = AR /40 sind tany , 


5 5 
' =¢'' = \R /4n sinw ad 3 
6 6 2 
AR, tan” (v + a4) 
"8° 8nsinv 





— aul cosy TT 
yt oe os - 1 ee ; 
7 X sin a, Z 


Sle Cie 


m-agpsm (cont. ) 
47Rj fi [or - ¥] Tt 
“4° y \osingy 4 
: ; 47R, fos [+ + ¥] au 
ne oe sin a 4 
t 
2 + 
o,. a T Oy 
g = @ +7 
5 8 : 
4n(R>-R,) sind 
i a a 
and 
3" - Aan (R> -R}) siny 
6 8 X 
For Y=rT 
=o_ + 
c=¢C o to) ,+4 T75, 
+2 _ 
TIF 4 cos ($, > 0 
where 
o a : 
7 t 
_ 2 2 
T >= TR, tan a, 
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re 
1! 

aN 
= 
a eas 
ae 
Yn) ~ 
5S 
le hr 
so 
i] 
MZ 
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m|a 





an Rk, ze 4nR5 
1) = -'—= =; and ) ON ee 
9 tan a, 2 I0 = Xtan a 


The formulas as listed above are used in Method II. To apply Method 
iit is necessary to delete all of the "cosine terms" from the above for- 
mulas (i.e. , each factor of cos (>; - ¢;) is a measure of the relative phase 
and should be replaced by zero). 


B.2.2 Computations for Configurations C, D, E, F, G, and H 
es ee ee 


Configurations C and D, being ogival in the front, were computed 
using the ogive formulas of Reference 5, which are listed below (b = 
maximum radius = 0.50 m and a = 1/2 nose angle)!: 


ot + cos ; 
el? Caeser oll 
sinza 
1 + } + cosa 
= (b cos a) 
sin a 


For 0< W<s -a 


¢ tan? a 
Co = 


1Smcos°y (1 - tan2a tan2 »)> 


‘An extensive discussion of the determination of the cross- section of an 
Ogive and a summary of Ogive experimental data will be contained in the 
next report on this contract. 
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4m tan (a/2) 


F z a< rae 
ors US 


ape 2 
ane P siny 





For Configurations E, F, G, and H, the techniques for C and D were 
employed in the vicinity of ¥ = 0 degrees, and the basic procedure usedon 
Configurations A and B (Method I) were employed in the vicinity of 
W =7/2-a,. The cross-section curve is faired in between the two points 
determined. 


B.2.3 Cross-Section Data at 1000 and 3000 Mc 


The data for Configurations I— VI, A, and B are given in Figures 
A.3-1—A.3-9. The cross-sections of Configurations A —H at 1000 and 
3000 Mc are shown in Figures B.2-2 and B.2-3. In examining Figures 
B,2-2 and B.2-3 it should be recalled that the cross-section of Config- 
uration I is approximately 107! mé throughout the aspect interval. 


B.3 CROSS-SECTIONS FOR CONFIGURATIONS A —H IN THE 200 TO 
ee a ee ee ee oe a 
500 Mc RANGE 


B.3.1 Configurations A -—F 


In Appendix A the nose-on cross-sections of Configurations I — VI 
were computed by a ''creeping-wave'' method. Applying this technique 
for Configurations A —F, the results shown in Table B.3-1 are obtained. 
The contribution from the rear, o2, is the same for Configurations 
A —F as it was for Configurations I — VI, but the nose contribution, ae 
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is in most cases less for Configurations A—F.' The determination of 

o, and o2 for each configuration yields, as it did in Appendix A for 
Configurations I— VI, an average cross-section, o] +02, and bound 
estimates, (Yo, tYo2)*. These three values are given in Table B.3-1. 
A comparison of these results with those for Configurations I — VI appears 
in Figures B.3-1 —B.3-3. 


Examination of these figures indicates that little can be gained, 
cross-section-wise, at 200, 350, or 500 Mc, by changing the re-entry 
shapes to those of Configurations A — F.’ 


B.3.2 Configurations G and H 


Although the cross-sections of Configurations A —F are not signifi- 
cantly less than those of the original Rudolph design, Configurations G 
and H, due to their ogival shape, can be expected to yield large reduc- 
tions in cross-section even at 200 Mc. 


B.3.2.1 Experimental Ogive Data at 200 Mc 


Experimental ogive data obtained by Ohio State University, and 
reported in Reference 15, included measurements made at 3000 Mc 


‘The reader is referred to Section A.2.1 for a discussion of this 
creeping wave method. Its application to Configurations A —F is similar 
to its use on Configurations I — VI in that the contribution from the rear, 
o2, iS approximated by that from a sphere of radius equal to 0.5 m for 
each of the Configurations A—F. The differences revolve around the 
determination of the contribution from the nose, o,. For Configurations 
I — VI the nose is spherical and thus go, in those cases is given by r Rn“ 
where R, is the radius of the spherical nose; for Configurations A -F 
the nose is pointed and thus the contribution, 0}, 18 given by 


Ne tan* a 
lor 


*This conclusion is also verified by the experimental data on the 
Carrot, the Cone, and the 7-OC models discussed in Section A.4. 
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(. = 10 cm) on two ogives: one of length 10 cm and maximum diameter 
2.05 cm, and the other of length 20 cm and maximum diameter 5.18 cm. 
Upon scaling this data to 200 Mc, one obtains ogives comparable in size 
and shape to Configuration G. Table B.3-2 contains a comparison between 
the dimensions of Configuration I (the unperturbed Rudolph shape), Con- 
figuration G, and these two ogives, denoted as "Ogive A" and "Ogive B". 


The data given in Reference 15 yield values of cross-section for 
the 10 cm ogive at four values of ¥(0, 34, 75, and 90 degrees), and for 
the 20 cm ogive at eight values of ¥(0, 15, 20, 43, 45, 50, 85, and 90 
degrees). 


Figure B.3-4 contains these data scaled to 200 Mc. In drawing these 
curves only the points defined by the angles mentioned above were used. 
Figure B.3-4 also contains the theoretical (physical-optics) curves for 


Ogives A and B. The theoretical curves donot differ from the exper- 
imental curves by more than a factor of 10 over most of the interval. 
: 


TABLE B.3-2 


COMPARISON OF DIMENSIONS 


(L = total length; D = maximum diameter; anda = 1/2 nose angle) 
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Ogive B: a@=29° L=3.00 m 
@—#@ Experimental (Horizontal Polarization) 
O——0O Experimental (Vertical Polarization) 
—— Theoretical (Physical Optics) 
Ogive A: @=23° L=1.50m 

@--@ Experimental (Horizontal Polarization) 
O-~-0 Experimental (Vertical Polarization) 
-~-- Theoretical (Physical Optics) 
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B.3.2.2 Cross-Sections of Configurations G and H at 200 Mc 


In view of the agreement between experimental and theoretical 
results for ogives shown in Section B.3.2, physical optics ogive methods 
were used for Configurations G and H at 200 Mc.! The results obtained 
are shown in Figure B.3-5, where they are compared with the Config- 
uration I data. 


B.4 CONCLUSION 


The perturbation analysis described above has indicated that, of the 
perturbations considered, the one referred to as Configuration G yields 
the largest reductions in cross-section in the aspect interval 0°< » < 60°. 
Capping the nose with a cone without adding any additional protuberances 
to the Rudolph shape does not decrease the nose-on cross-section of the 
warhead to any appreciable extent. Capping the nose with a cone and 
rounding-off the step results in large reductions in the noSe-on cross- 
section at the larger frequencies, but at 200 Mc no appreciable reduction 
is obtained. If, in addition to capping the nose of the warhead with a 
cone and rounding-off the step, the rear of the warhead is also capped 
with a cone (as in Configuration G), then comparatively large reductions 
in the nose-on cross-section can be obtained at 200 Mc also. This is 
illustrated in Figure B.4-1 for Configurations A, C, E, and G, compared 
with Configuration I. The other four perturbations (B, D, F, and H) are 
not included in this graphic comparison, since they would not enclose 
Configuration I. 


It is seen from Figure B.4-1 that the nose-on cross-section of Con- 
figuration G is smaller than that of Configuration I by factors of approxi- 
mately 10° at 3000 Mc, 107 at 1000 Mc, and 102 at 200 Mc. 


'The use of these ogive methods at small wavelengths has been 
shown to be valid by the work of Sletten (Ref. 14) and by the work pre- 
sented in Reference 40. Since there is some contradictory evidence in 
recent O.5.U. work, additional discussion on this point will be made in 
the next report under this contract. 
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APPENDIX C 


CAMOUFLAGE MATERIALS 


A radar absorbing material has been developed by the Deutsche 
Magnesit,AG,Munich, Germany. Representatives of this company were 
invited by the Rome Air Development Center to fabricate their absorbing 
material at RADC and then apply it to models furnished by RADC. The 
results are very good when the material is manufactured for use in one 
small frequency region lying in the S- to K-band region. In fact, it was 
claimed that there was a 26 db reduction in cross-section at X-band. 


The German scientists had no mass production methods available for 
the fabrication of their material. They also had no method available for 
making broadband materials, but they thought they could do so without 
too much difficulty. They attempted to do this while at RADC. The re- 
sults were not good, but this may be because of the insufficient time 
available. 


After they had left, work was started at RADC under the direction of 
J. Vogelman, Director of Electronic Warfare Laboratory, to find out 
why the obtained reduction existed and to analyze its division into absorp- 
tion and increased scattering at other than the back-scattering direction. 
It was found that half the decrease in db was due to absorption and the 
other half was due to off-angle scattering. 


Due to the newness of this material, and due to the fact that broad- 
band techniques have not yet been completed, the only conclusion to be 
drawn at this time is that the Deutsche Magnesit material may be usable 
for camouflaging missiles and aircraft. Since RADC is continuing its 
investigation of this material for the U. S. Air Force, it is felt that its 
analyses will show how good or bad this material is. Atthis time the 
authors do not recommend any further action than that which is presently 
going on at RADC. It is recommended that people in the ICBM field fol- 
low RADC's latest results on this material for possible applications to 
ICBM. 
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APPENDIX D 


UTILIZATION OF FISSION WASTE PRODUCTS FOR DECOY PURPOSES 


D. 1 INTRODUCTION 


Fission waste products are radioactive and therefore produce ion- 
ization. It has been conjectured that if the wastes are carried in a decoy, 
the radar cross-section of the decoy might be sufficiently enhanced to 
improve the chance of the decoy being confused with the primary vehicle 
by proximity-fused interceptors. Also, if the decoys were used to draw 
off enemy defenses to a false location by confusing ground-based radars, 
one would want to increase their detectability. This idea is considered 
here for two types of missiles and found impracticable. Another group 
has reached the same conclusion (Ref. 16) for sowing the wastes from 
airplanes as a sort of chaff. 


The first type of missile considered is a hypothetical intercontinental 
ballistic missile. The following properties, which should be duplicated 
by an effective decoy, are assumed: a velocity of 10 km/sec, and flight 
altitudes ranging from ground level to above 200 km. Since in its nose- 
on descent such a missile would not have a large cross-section, it is 
assumed that something between rm and 10 m square meters would be a 
reasonable radar cross-section to simulate. 


The second type is a V-2, representative of the shorter-range mis- 
siles. Experimental data quoted in Reference i7 indicate typical altitudes 
up to about 90 km and velocities ranging from about 3.25 km/sec down 
to about 0.5 km/sec at different altitudes. The size of the cross-section 
to be simulated is the same as above. 

D. 2 PHYSICAL PRINCIPLES INVOLVED 
D.2.1 Critical Density 


Upon solving the propagation equations in an ionized medium for a 
given radar frequency, it is found that if the electron density exceeds a 
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critical value the surface bounding this dense region acts as a perfect 
reflector (Ref. 18). This criticaldensity is given(Ref. 1)by Ne = (f2/81) x 
10electrons/cm3, where f is the frequency in megacycles. A sub- 
critical-density region surrounding the critical-density contour acts as 

a diffuse reflector, reducing the back scattering from the critical-density 
region. (The cross-section of a region which nowhere attains critical 
density is very much less than that of a geometrically identical critical- 
density region.) The critical-density region should be at least a few 
wavelengths thick if the propagation equations are to have meaning. A 
radius of about one meter should be sufficient, even down to L-band. 


D.2.2 The Nature of the Fission Waste Products 


The radioactivity of the fission waste products consists almost ex- 
clusively of B- and ?-radiation. Therefore, other types of radiation do 
not come into consideration. A representative figure for the power re- 
leased is 100 watts/gm, divided equally between the 6's and Y's. This 
figure, based on radiochemical analysis, does not take into account the 
loss of output due to self-absorption in a bulk source. The mean energy 
of both types of radiation is 0.7 Mev, the maximum B-energy about 5 
Mev (Ref. 191). 


D.2.3 Range-Energy Relations For £'s and 7's 
I Pe ee 


For electron energies such that v>c (since v/c = 0.5 corresponds to 
less than 80 KV, this covers nearly all of the B-activity present here), 
the energy loss of a 6 per g cm™~* of absorber is given by (Ref. 20) 





3/2 
TE 3/ 


2 2 
k = 4 N(Z/A) ar.* mz c*| In —\ -1,45 
_ * BE, 1/2 yz) 


? 


1The authors wish to express their thanks to J. V. Nehemias for 
communicating these results to them. 
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where N = Avogadro's number 


atomic number 


Z/A = atomic weight ° for the stopping element, 
Ye = e*/me cf, the classical electron radius, 


electron rest mass, 


5 


total electron energy, 


ie 
J 
r 


2 


E.. = rest energy, m,c”, and 


I(Z) = ionization potential of stopping material. 


Since k varies monotonically by less than 30 per cent from 0 up to 5 Mev, 
it may be taken as constant at its value for the mean 6 energy of 0.7 Mev 
with little error. Using as order of magnitude figures for the fission 
products Z/A = 0.45 and I(Z) = 8 ev gives k = 1.6 Mev/g cm~-, Using a 
density p = 5g cm ~ for the wastes yields a loss rate of 8 Mev/cm, ora 
range of about 1 mm for 0.7 Mev 8's in the source. For ranges in air, 
Z/A = 0.5 may be used, and little error is involved in taking I(Z) = 13.5 
ev, correct for atomic oxygen. p must correspond to the altitude con- 
sidered. 


The range of the electrons in the waste material being so small, ab- 
sorption of the electrons in the source itself (self-absorption) will greatly 
reduce the power output of the source into the air, allowing only those 
electrons originating near the surface to escape. For electrons of given 
energy, consider a strip of area A and thickness t equal to the range, and 
compute what fraction of their energy will escape through the outer side 
of the strip (averaging over points of origin). This fraction is F=1- (s/t), 
where s is the mean path of the electron in the strip (since energy loss 
is proportional to path length). 
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cos~!(a/t) 


Ce 


F =—————_— Ada singde@(t-a/cos 6) 


t 

] l ] 

ae | da [tt -a/t)+aln 2/1) cake - t¢/2 2/4 | = a 
0 


This derivation involves the assumption that the thickness of the strip is 
small compared to its radius of curvature (which agrees with the actual 
dimensions considered later), and assumes a straight-line path. Ifa 
thicker strip is used, only electrons originating within a distance t from 
the surface will contribute. Neglecting range straggling, elec- 

trons originating at a distance from the surface greater than their range 
will be absorbed before they reach the surface. The use of a thinner 
strip is not advisable, since the reduction of surface activity due to the 
presence of less source more than counterbalances the lowered self- 
absorption. 
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For a 0.7 Mev electron in air, k = 1.5 Mev/g cm™*. Near sea level, 


the air density is about 1073 g/ cm?, Thus the energy loss per unit path 
length is 1.5 x 10-3 Mev/cm, and the range is about 500 cmor5m. As 
the altitude increases, the air density drops off more and more rapidly 
(e.g., 10712 g/cm? at 170 km) and the range increases proportionately. 


Gamma-rays can be considered as attenuated as exp(-px) (at least, 
for px of order unity or less). For the ¥-energies of interest, the inter- 
action of the gamma-rays with matter is practically all due to the Compton 
effect (pair production competes prominently only at much higher energies, 
the photoelectric effect only at very much lower ones). The absorption 
coefficient, p», for the Compton effect is (Ref. 21) 


Zo 2f l+al 2c 
pe 2m Np i re { 43] Ha) - bn +20)] 





] 1+ 3a 
] me 
+ 35 In (1+ 2a) ep 


Where a = ¥ -energy in units of the electron rest mass. 


Bn 


3 


For a 0.7 Mev y-ray in air near sea level, p = 4x107? cm 
hence, the Y-ray loses~ 63 per cent of its energy in about 300 meters. This 
distance varies with air density (consequently with altitude) in the same 
manner as the electron range. Thus, the ¥-rays are about 100 times 
less efficient than the 6's in producing ionization. 


Inasmuch as the available power is evenly divided between B's and 
y's, and that portion of it in the form of 6's is far more effective, con- 
sideration is limited to the electrons alone. ! 


In view of the high self-absorption of the source for electrons, the 
power output per gram quoted above (neglecting absorption) is not di- 
rectly applicable to a bulk source. Instead, the useful quantity is the 
power output per unit area of the source. There are 50 watts/gm emitted 
in the form of B's. Using a density p=5 g/cm? for the waste products, 


'The overly generous assumption that all the radiation is in the form 
of B's would reduce the source size by less than a factor of two, leaving 
the conclusions unchanged. 
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as before, this is 250 watts/cm?, From the discussion of self-absorption 
above, the power released by a layer of thickness t cm, equal to the range, 
is 1/8 X 250 Xt watts/cm*, and deeper-lying material does not contribute. 
Since the rate of energy loss is 8 Mev/cm, t = E/(8 Mev/cm) and the out- 
put is (250/64) E watts/cm“, where E = electron kinetic energy in Mev. 
Averaging E over the electron spectrum gives 2.7 watts/cm%. This is 

the surface activity to be expected from a thickness equal to at least the 
range of the most energetic 8's (about 6mm), A thinner strip will give 
less surface activity and a thicker one will yield no more. 


D. 3 COMPUTATIONS 
D.3.1 Introduction 


This section contains a computation of the source strength required to 
maintain a critical density of ionization over a volume adequate for decoy 
purposes. It will be shown that an impractically large source is needed. 
Since this consideration alone rules out the use of fission wastes, others 
need not be treated quantitatively. However, these additional limitations 
are stated here. 


First, there is the fact that the B-energy spectrum is continuous. 
This implies a continuous distribution of ranges, and therefore an ion 
density distribution which tails off away from the source. Thus, beyond 
the region of critical density a subcritical density is present which re- 
duces the radar cross-section. 


Further, the phenomena of recombination and attachment to form 
negative ions have been ignored. These would work to lower the electron 
density and require a larger source. On the other hand, if recom- 
bination is not sufficiently rapid behind the moving source, an ion 
trail may be formed. In the latter case, the radar reflection from the 
ionization produced by the decoy would not resemble a moving mis- 
sile. 


The fact that the B's must be relied on to produce the ionization 


raises a problem as to mounting the source on the decoy. As was in- 
dicated in the discussion of self-absorption, one centimeter of solid 
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material would stop all the 6's. This effectively precludes carrying the 
source inside the decoy. 


The nature of the fission waste products raises serious handling 
problems. The high level of radioactivity would make the loading of the 
decoy unpleasant and result in contamination difficulties at the launching 
site. The effects of this radioactivity upon the chemical fuels powering 
the decoy should also be considered. Finally, the large amount of power 
dissipated in self-absorption would create significant thermal problems. 


D.3.2 Source Strength Requirements for Critical Density 


To simulate a missile, ionization of at least critical density must be 
made to fill a volume comparable with the volume of the missile. Hence, 
at any moment, a critical density region of radius r, (ideally a couple of 
meters) is required. As the decoy moves, the region that has to be filled 
changes. If the missile velocity is L cm/sec, a cylinder. of radius r) and 
length L cm must be filled in one second. Since the source is located on 
the axis of the cylinder, the ion density decreases as the distance from the 
axis increases. Therefore, critical density must be achieved near the 
surface of the cylinder (which will result in a higher density being pre- 
sent further inside). For simplicity, assume the source to be a sphere 
of radius rg (this shape also minimizes reabsorption, in one part of the 
source,’ of electrons that have escaped the surface at another point). To 
fool the radar, the ionized layer must be at least a meter thick, requiring 
that r) - ro #100 cm. 





Given a source whose power output per unit surface area is 2.7 
watts/cm4, subject to the conditions just enumerated, the minimum ro 
required will now be computed. The calculation is not carried out ex- 
actly; instead, a lower bound to the answer is obtained in three different 
ways (which is the most restrictive depends on the altitude and velocity 
considered). 


D. 3.2.1 ICBM (near sea level) 


To maintain the critical density out to rj, it is necessary to supply, 
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every second, a cylindrical shell of radius r; and thickness dr with an 
amount of energy 


€=2nrr,; Ldr DI(Z) 
where D = critical density, and 


I(Z) 
Let f(E) dE 


ionization potential (13.5 volts). 


number of electrons supplied by the source, with 
energy lying between E and E + dE, 


k(E) = energy loss of electrons per unit path length (approx- 


imately independent of BE), 


ds 


path length of electron through the cylindrical shell 
(ds 2 dr), and 


E' = minimum energy of electron at source if it is to 
have a range sufficient to arrive at the cylindrical 


Shell. 


Then 


° ©) 
é=k { f(E)dE ds 
E! 


CO 


= « | BIE = 
E 


FE! 


& 
ds 
ak == Ef(E)dE 
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gas, | EME)dE. 
0 


For different orientations of the electron path, ds is proportional to the 
path length traveled in order to reach the cylindrical shell (minimum 
range), which in turn is proportional to the minimum required energy E'. 
Hence, ds/E' is independent of direction and can be replaced by dr/E), 
where E; is the minimum energy for radial motion. If the dimensions 
of the source are small compared with r,, E, = kr). Hence, 


ao} 
é€<,% 8 | Ef (E)\dE 
E! 
0 


— {(E)d 
E] Ef(E)dE 
0 


d 
— Ef(E)dE 
ry 


le, 8) 
and | Ef(E)dE > 2nr 
0 


1° LD I(Z). 
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o,@) 


The energy supplied, | Ef(E)dE, comes from a sphere of radius ro. 
0 


In one second, this source supplies 2.7 joules per cm® of area or a total 
of 2. 7X 4nr )4 joules = 2.12 X 1029 Gr ev. For minimum energy input for 
the ICBM (L = 10° cm), 


2.12 x 1029 rp? ev = 2mr,2 10°cm x 6.9 x10!! x 13.5 ev 


ro¢ 
——— = 0.28 
ry? 

Ty = 0,53 ry: 


Since r, - Xo 7 100cm for simulation, the minimum rg = 110 cm, so the 
diameter of the sphere has to be more than two meters. Thus at low 
altitudes, the fission products ionization source is at best no more ef- 
fective as a decoy than a metallic object of comparable size would be. 


Neglect of the source dimensions is not essential. 
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The path difference of an electron arriving at a point on the cylinder 
from different places on the source can be taken into account. Assuming 
isotropy of emission direction from the source, the correction factor is 


max 
sin8d6@ r 


far oo Pat Renee 
Tijav. F] 


max 
sin@d@ 
0 
ay an 
From the geometry, r“-r *+r4%-2r 7. cos @ and @ = cos”! , 
] 0 0 1 max ry 
r 
Let x =—, cos 9 =u, then 
r 
] 
— =Vlt¢xt-2xu 
= 


] 

VA 
| Vi+x - 2x udu 
x 





1/2 
_ (1 ae fee : a-x3/?] 


For x = 3 


Gi = 0,7, 
| 
av 


The minimum diameter of the sphere is changed to 1 1/2 meters. The 
conclusions are unaltered. 
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D. 3.2.2 ICBM (higher altitudes). 


At low air densities where the range of most of the source electrons 
greatly exceeds r), the replacement of ds/E by ds/E' represents a 
gross underestimate of the source strength required. Instead, 1/E is 
left inside the integral and a mean value of ds, or better ds/dr, is used: 


ds S ] 


dr- ry cos 9 





Electrons will leave the cylinder at angle 0 ranging from zero toa 
maximum value such that s is the electron range. An overestimate is 
obtained by choosing s,,a, = Eas) k and assuming an equal number of 
electrons leaving at all values of 8 in the range so defined (i.e. , 0<€ O0< 


Omax = cos71 (kr / Eas) ). 





Then a 40 
ds . cos®@ 4 | 
=e Onar COS Gs as n (tan Q@may + sec 9max) 
i. 
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where E is the mean energy. 
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and, at the minimum, 


2 1.89 Xlo-8 Lr, 











r = 
: 6. 25 
Pin 
pry 
] 
The condition kr; « E,,,, means 1. 6 x 106 pr); <«< 5x 10° or ry >). 
Z5 
Thus, ln =ln DS gee a ay ee and 

pry p 6.25 p 


r2 1,89x107e8 7 


Yl 7 o inl 
p 


Near sea level, p= 107? em/cm? and a reasonable r, is 200 cm, so 
that kr}/Emax = 0.06. At higher altitudes, the inequality is even better 
satisfied. For the minimum rg, impose again the condition that 

r; - ro = 100 cm (minimum ionization layer for simulation), 


ry? 1.89 x10-8 LZ 
a Ae aoe 
ry t 100 eee 


For the ICBM at different altitudes, the source size is given below: 


Altitude in km Air Density in gm/cm> (Ref. 22) 


4.21xX 1074 0.5m 
















9.34% 1079 


1.16 x 107° 
-10 


1.2m 





26 m 







8.60 x 10 21 km 


7.5104 km 
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Thus, while this lower bound on the source is somewhat less strin- 
gent than the previous one at low altitudes, it imposes far greater (in 
fact, quite unattainable) limitations at the higher altitudes. 


It would appear at first sight that the large value of rg obtained at 
the higher altitudes conflicts with the omission of source dimension 
considerations from the derivation. This is not correct. Indeed, in- 
cluding these considerations results in changing the maximum of s to 
(Enax/kK) + Pq. This simply makes 


r 
Emax +p, Emax (1+*%o __ 
k Emax 


Since krg< kr, «<E , this correction factor can be neglected. 


Max 


D. 3.2.3 The V-2 


The V-2 has a lower velocity than an ICBM. Furthermore, its 
maximum altitude is only about 150 km, which is considerably less than 
that of an ICBM. These differences reduce somewhat the power require- 
ments. Applying the calculation of Section D. 2, 


Altitude in km 48 63 88 


Velocity in km/sec] 0.995 1, 32 1.08 


3.25 
(Ref.17) 


9 


Density in gm/cm |] 9.341079 1.48 x10°°| 2.42x10! 
(Ref. 22) 


5.65 X10. 


Diameter 16 m 1.1 km 
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As the altitude increases, ®y9 x reaches a cutoff value due to the 
size of the cylinder, namely Ona, = tan7! (L/2r)). 


L/2 





Then q 1 
dr/ay max 


] L L \2 
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For (L/2 2 > rane this leads to a maximum value: 
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2in & 
ron | .9 x 10-8 
ry 9 


Set r) = ry + 100 cm (minimum): 


2 Cr 
ro” In rg + 100 1.9X1078 7 
ry + 100 p 


Results obtained in this way are: 


Altitude in km 


Diameter in m 





At an altitude of 88 km, this approach becomes meaningless because 
ro > L. 
D.4 CONCLUSIONS 


The proposal to utilize fission waste products in decoys to simulate 
the radar cross-section of conventional missiles or the ICBM has been 
examined. From source strength (and size) considerations alone, the 
idea has been shown to be impractical. Even near sea level, a lower 
limit on the power required leads to a source of size comparable to that 
of a metallic object giving the same cross-section. At higher altitudes 
(over most of the range for the V-2 and practically all of the range for 
the ICBM), the source required is unattainably large. Other difficulties, 
mentioned but not treated quantitatively, raise further doubts as to the 
feasibility of the proposal. 
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APPENDIX E 


DECOYS 


Before considering possible defenses against decoys from the radar 
cross-section point of view, we will summarize the work which has been 
done for the aircraft case. 


One first computes the radar cross-section of the aircraft. Then 
one investigates means of simulating its radar cross-section and flight 
characteristics with a smaller body. In Reference 5, the radar cross- 
sections of the B-47 and B-52 were computed from 65 Mc through X-band. 
It was found that an object of approximately 1/100th the volume of the 
aircraft could be used to reproduce the radar cross-section; the object 
could also reproduce the flight trajectory for one flight. In Reference 2 
the possible defensive methods outlined below and analyses for the design 
of decoys to simulate the reflection characteristics of the aircraft were 
presented. In Reference 6, the scattering matrix of the B-47 was deter- 
mined at S-band. The latter was useful in discussing polarization meth- 
ods of distinguishing between decoys and aircraft. If one assumes that 
the defense has linearly polarized monostatic radars of fixed frequency 
(L-band or higher), then the radar cross-section of the aircraft can be 
simulated through the.use of corner.reflectors in decoys. 


However, the enemy may then use frequency comparison methods to 
determine whether the vehicle is an aircraft or a decoy. For example, 
if the enemy's main radar is at S-band and its frequency comparator 
radar is at X-band, one would expect the cross-section of the aircraft 
to remain approximately the same for the two radars while the cross- 
section of the decoy would increase by a factor of 10 from S-band to 
X-band, due to the fact that the cross-section of a corner reflector de- 
pends on wavelength approximately as 1/r2 . On the other hand, it has 
been our experience (Ref. 23) that the radar cross-sections of aircraft 
such as the B-47 and the F-86 are only slowly varying functions of 
wavelength as soon as the wavelength becomes small with respect to the 
major dimensions of the aircraft. 
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This detection scheme can be countered as follows. Corner reflec- 
tors with appropriately curved faces will simulate the aircraft at S-band 
but cannot be recognized by the above technique. This is so because 
curvature of the faces of a corner reflector greatly modifies the behavior 
of its cross-section as a function of \. A feeling for this can be gotten 
by observing that the angle error formulas for corner reflectors are 
wavelength-dependent, and thus changes in angle which are insignificant 
at one wavelength can cause severe reduction in the radar cross-section 
at smaller wavelengths. It has been shown that a corner reflector with 
spherical faces can be designed that eliminates the 1/\4 wavelength de- 
pendence. The beauty of the spherical face idea is that after one matches 
the cross-section at one frequency it remains essentially a constant for 
all higher frequencies. 


Another means of ''unmasking" corner reflectors in decoys is the use 
of a transmitter which is circularly polarized (e.g., to the right), to- 
gether with two receiving antennas, one of which is right-circularly po- 
larized and the other of which is left-circularly polarized. When the cor- 
ner reflectors are so oriented that single or triple reflections dominate, 
the left-circular polarization receiver will get 10 to 20 times as much 
power as the right polarized receiver. When double reflections dominate, 
the reverse Situation applies. In any practical case, one of these orien- 
tations will occur. Thus, decoys using corner reflectors could easily 
be differentiated from aircraft (for which the power distribution between 
the two receivers varies between 3:2 and 2:3). This method of detection 
can be countered by the use of dielectric material on one of the faces of 
the corner reflector (Ref. 24). The procedure consists of properly coat- 
ing one of the faces with a dielectric, so that the ray picture remains 
essentially unchanged, but the polarization is rotated in such a manner 
that the return is now distributed between the two receivers in roughly 
the ratio expected from aircraft. Another means of defeating this method 
of detection is by the use of holes of appropriate size in one of the corner- 
reflector faces or the use of one-half-wavelength protuberances. These 
last two approaches are not as efficient as the Ohio State University tech- 


nique (Ref. 24) because they cause a significant degradation in the energy 
returned. 
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Another possible means of discriminating between decoys and air- 
craft would be the use of bistatic radars (Ref. 25 and 26). Bistatic radars 
are presently being used in the McGill Fence and the DEW line. Dis- 
crimination is based, in this case, on having one receiver near the trans- 
mitter and one or more receivers at a distance from the transmitter, In 
the case of aircraft (assumed to fly between the remote receiver and the 
transmitter), the radar cross-section at the remote site is often much 
larger than, and in almost all cases at least comparable to, the radar 
cross-section at. the receiver which is near the transmitter. In the case 
of decoys, the cross-section at the transmitter site is always much 
larger than the cross-section at the distant receiver, due to the property 
of the corner reflector of reradiating energy into the quadrant from which 
itcame. Some tentative schemes (obviously not fully adequate) have been 
discovered to counter this method of detection (Ref. 2). 


Two other methods for the detection of decoys have been brought 
forward; one involves scintillation and glint and the other broadside dis- 
crimination (it is assumed in the discussion that the frequency of de- 
tection is above 1000 Mc). The power spectra for one-third-scale V-2 
type ballistic missiles are contained in this report. The authors have 
no reason to believe that the power spectra which would be obtained for 
the Rudolph missile would be significantly different from these. In fact, 
if one analyzes the return from aircraft like the B-52 and B-47 (Ref. 27), 
he finds the shapes of the power spectra are rather similar to those for 
the missiles. This, of course, does not mean that special equipment can- 
not be built to detect spectral differences. However, it would be extremely 
difficult to do this without a thorough knowledge of the structural behavior 
of the enemy's missile and the enemy's decoy. The broadside discrim- 
ination method is more applicable to the aircraft problem than the mis- 
sile problem, since in the search phase the Rudolph missile would be out 
of the sensible atmosphere, and as a result Rudolph never would have to 
"put its broadside forward.'' In fact, one of the best methods of reduc- 
ing the cross-section of the Rudolph missile would be to guarantee that 
its nose is always pointing in the general direction of the enemy's search 
radars, with the temporary alteration of the metallic skin described in 
this report applied, so that the cross-section of the Rudolph warhead in 
the nose-on direction would be significantly reduced. 
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In the case of low frequency detection (e.g., the 60 to 100 Mc 
range), the decoys must carry barrage jammers if they are to duplicate 
the low frequency cross-sections of the Rudolph missile. However, the 
background due to cosmic noise, meteors, and stray ionization sharply 
limits the feasibility of detection by search radars in that range. It has 
been said that if American radars were operated at 60 Mc, it might be 
possible for the "Commissar of Jamming of the USSR" to push a button 
and create so much radiation that he could jam out all our receivers in 
this frequency range. Due to the great range of propagation at these 
low frequencies, it is clear that the above statement has an aura of truth 
about it. 


The bistatic method for detecting decoys is not presently being count- 
ered. However, the data-gathering and data-processing equipment nec- 
essary for the use of this method would tax our national economy. Al- 
though it has been assumed that the United States has more computing 
equipment and data-processing equipment available than the Russians, it 
is felt that even the United States could not afford in the next few years 
to use a great many remote receivers and the associated data-process- 
ing equipment in order to discriminate between decoys and aircraft. 

Thus, it seems reasonable to assume that the USSR also will not be able 
to afford such defensive procedures. It appears, therefore, that decoys 
could substantially increase the "probability of kill" of the ICBM. Of 
course, until one is sure that the enemy can defend against such missiles 
the necessity for a decoy effort would not seem to be established. How- 
ever, if one waits until the knowledge is obtained that a defense is pos- 
sible, it might easily be too late to design and build the necessary number 
of decoys. 


3 


Decoys, of course, can be used for both the long-range and short- 
range efforts; they can be used against search radars and against local 
defense radars. It may be possible to store decoys in the tanks section 
(if it is not far removed from the warhead) which can be used against 
elther set of radars. It is possible to make decoys which are aerody- 
namically more highly performing than the Rudolph warhead (they 
might even catch up with and pass the warhead in flight after re-entry 
into the atmosphere). 
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APPENDIX F 


FURTHER ANALYSIS OF THE UNIVERSITY OF MICHIGAN 
DROP-TEST DATA (POWER SPECTRA FOR EXPERIMENTAL DATA) 


F.1 INTRODUCTION AND SUMMARY 


As stated in Section IV and in Appendix E, one method the defense 
might attempt to employ to discriminate between decoys and aircraft is 
based upon possible differences in the fine structure of the radar return. 
An investigation of this possibility involves a thorough study of the fine 
structure of the radar return from aircraft, decoys, and missiles. 
Others have expressed the belief that power spectra associated with the 
radar return from these vehicles would be of significance in this detec- 
tion problem. 


Power spectra for manned aircraft have previously been studied 
(Ref. 27); to the best knowledge of the authors, the power spectra pre- 
sented here are the first that have been computed for the radar return 
from missiles. The experimental data involved in the present analysis 
are the measurements made during The University of Michigan's drop 
tests reported in References 1 and 3. These computations were spon- 
sored by the Bendix Aviation Corporation under Purchase Order No. 
S-96527. 


The one-third-scale V-2-type models involved in the experiments 
are referred to here as Objects 59, 61, 62, 64, and 65, The first three 
models had swept-back fins and the last two had rectangular fins. 


Three separate items were determined under this study: a frequency 
distribution of the radar cross-sections of each model at each wavelength; 
the autocorrelation function for each run; and the power spectrum for 
each run. (The procedures employed and the data obtained appear in 
Sec. F.2—F.4.) 


The study brings into focus two basic questions: first, whether or 


not a fine structure analysis for the missile-decoy problem can be ex- 
pected to provide significant information; and second, the question of 
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whether or not the power Spectra are the correct quantities to consider 
in this analysis. This study deals with finite samples. The classical 
theory is based upon a sample of infinite size. A certain amount of arbi- 
trariness is involved when the theory is modified for application to a 
finite sample. Alternative definitions of the power spectrum and auto- 
correlation function have been given in the literature (Refs. 28 and 31). 
It is pointed out that these definitions are definitely not equivalent (Sec. 
F, 5). 


The power spectra obtained by the procedure outlined in this appen- 
dix are displayed graphically in Section F.4. Even though the main con- 
tribution to the radar cross-sections of V-2-type missiles usually comes 
from the fins, the variation in the shapes of the curves found for Objects 
59, 61, and 62 (swept-back fins) appears to be as great as the variation 
between these and Objects 64 and 65 (rectangular fins)’, 


It is possible that a study of a larger sample of missile runs by the 
approach used here might begin to display a pattern of differences be- 
tween power spectra. An attempt to obtain greater discrimination by 
lengthening the time of observation for a given run would have no bearing 
on the practical defense problem, since the times available are even less 
than the intervals of observation used in the drop tests. The practical 
problem is further restricted, if a pulse-modulated radar is used, since 
the accessible frequency interval is limited; for the drop-test data, the 
upper limit was approximately 400 cps. 


some other approach to the question of how the infinite sample theory 
can be modified for application to a finite sample might conceivably yield 
Spectra which are more sensitive to differences between radar returns 
for different objects. 


If one is convinced that the fine structure of the target return is the 
most important information available for discriminating between aircraft 


‘In fact, there appears to be little difference between the spectra 
found here for missiles and those obtained for a B-47 aircraft (Ref. 27). 
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and decoys, it seems clear that much further analysis is necessary 
before one can hope to instrument this concept into a feasible discrimi- 
nator between decoys and aircraft. Such an investigation should include 
an analysis of the mathematical functional representation of the fine 
structure to yield the maximum amount of information. After the above 
method of representation is established and a method of presentation is 
fixed, one should determine (probably experimentally) whether the dif- 
ferences in the fine structure between missile and decoy are greater 
than the differences between missile and missile. The latter differences 
are expected because of differences in the conditions under which the 
observations are made or unanticipated alterations in enemy missiles 
(as to structure or structural behavior). If this investigation is to be 
used to analyze vehicles which are out of the atmosphere for almost all 
of their trajectories, one should bear in mind the possibility that tran- 
sient effects (which disappear rapidly in flights through only the sensible 
atmosphere) might not be damped out. It is clear that, for economic 
reasons, vehicles cannot be fired out of the atmosphere to analyze tran- 
sient effects which, by their very nature, may not be reproducible from 
air frame to air frame. 


In short, one should first complete a theoretical investigation to 
determine the best method of collecting fine structure data, assembling 
it, representing it, and then displaying it. When this is done, one should 
see whether, after considering the physics of the ICBM problem, the 
observables are really significant and physically meaningful in the opti- 
mal theory. 


F,2 THE FREQUENCY DISTRIBUTIONS 


The frequency distributions obtained for each of the ten runs appear 
in Tables F.2.1—F.2.10. Each value of cross-section was computed 
by averaging the power received over consecutive groups of six pulses 
each. The aspect angles for the measurements contained in these tables 
are such that the aspect angle never exceeds 15 degrees off-nose and in 
most cases is between 5 and 10 degrees. 
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TABLE F,2-1 


FREQUENCY DISTRIBUTION — OBJECT 59(S-band) 





o Interval No. of values 
(m“) in interval 

0.00 - 0.01 208 
0.01 - 0.02 107 
0.02 - 0.03 45 
0.03 - 0.04 20 
0.04 - 0,05 12 
0.05 - 0.06 8 
0.06 - 0.07 16 
0.07 - 0.08 26 
0.08 - 0.09 6 

448 
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TABLE F.2-2 


FREQUENCY DISTRIBUTION — OBJECT 59(X-band) 
a NG 


o Interval No. of values o¢ Interval No. of values 
(m“) in interval (m*) in interval 
0.00 - 0.01 45 0.17 - 0.18 4 
0.01 - 0.02 50 0.18 - 0.19 5 
0.02 - 0.03 55 0.19 - 0.20 5 
0.03 - 0.04 54 0.20 - 0.21 5 
0.04 - 0.05 48 0.21 - 0.22 4 
0.05 - 0.06 28 0.22 - 0.23 3 
0.06 - 0.07 28 0.23 - 0.24 4 
0.07 - 0.08 16 0.24 - 0.25 3 
0.08 - 0.09 19 0.25 - 0.26 ] 
0.09 - 0.10 19 0.26 - 0.27 ] 
0.10 -0.11 16 0.27 - 0.28 4 
0.11 -0.12 10 0.29 - 0.30 3 
0.12 - 0.13 5 0.31 - 0.32 l 
0.13 - 0.14 5 0.59 - 0.60 
0.14 - 0.15 5 0.64 - 0.65 ) 
0.15 - 0.16 9 "— 
0.16 - 0.17 I 


SECRET 


UNIVERSITY OF MICHIGAN 
2428-3-T 


TABLE F.2-3 


FREQUENCY DISTRIBUTION — OBJECT 61 (S-band) 


o Interval No. of values 

(m*) in interval 
0.00 - 0.01 40 
0.01 - 0.02 44 
0.02 - 0.03 27 
0.03 - 0.04 16 
0.04 - 0.05 8 
0.05 - 0.06 13 
0.06 - 0.07 14 
0.07 - 0.08 22 
0.08 - 0.09 17 
0.09 - 0.10 6 
0.10 - 0.11 5 
0.11 - 0.12 4 
0.12 - 0.13 3 
0.13 - 0.14 2 
221 
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TABLE F.2-4 


FREQUENCY DISTRIBUTION —OBJECT 61(X-band) 
————— ee OES NG] 


o Interval No. of values 


(m“) in interval 
0.00 - 0.01 49 
0.01 - 0.02 | 47 
0.02 - 0.03 17 
0.03 - 0.04 | 22 
0.04 - 0.05 12 
0.05 - 0.06 12 
0.06 - 0.07 5 
0.07 - 0.08 4 
0.08 - 0.09 3 
0.09 - 0.10 10 
0.10 - 0.11 9 
0.11 - 0.12 6 
0.12 - 0.13 8 
0.13 - 0.14 4 
0.14 - 0.15 5 
0.15 - 0.16 3 
0.16 - 0.17 l 
0.17 - 0.18 3 
0.18 - 0,19 2 
0.19 - 0.20 1 

223 
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TABLE F.2-5 


FREQUENCY DISTRIBUTION — OBJECT 62(S-band) 


o Interval No. of values 
(m“) in interval 
0.00 - 0.01 123 
0.01 - 0.02 107 
0.02 - 0.03 117 
0.03 - 0.04 62 
0.04 - 0.05 55 
0.05 - 0.06 48 
0.06 - 0.07 46 
0.07 - 0.08 21 
0.08 - 0.09 14 
0.09 - 0.10 12 
0.10 - 0.11 15 
0.11 - 0.12 3 
0.12 - 0.13 6 
0.13 - 0.14 4 
0.14 - 0.15 4 
0.16 - 0.17 2 
0.17 - 0.18 = 
641 
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TABLE F.2-6 


FREQUENCY DISTRIBUTION — OBJECT 62(X-band) 


o Interval No. of values o Interval No. of values 
(m“) in interval (m*) in interval 

0.00 - 0.01 152 0.15 - 0.16 3 
0.01 - 0.02 115 0.16 - 0,17 6 
0.02 - 0.03 80 0.17 - 0.18 2 
0.03 - 0.04 74, 0.18 - 0.19 3 
0.04 - 0.05 43 0.19 - 0.20 2 
0.05 - 0.06 35 0.20 - 0.21 5 
0.06 - 0.07 25 0.21 - 0.22 3 
0.07 - 0.08 13 0.22 - 0.23 Z 
0.08 - 0.09 15 0.25 - 0.26 ] 
0.09 - 0.10 19 0.27 - 0.28 ] 
0.10 - 0.11 12 0.35 - 0.36 ] 
0.11 - 0.12 8 0.36 - 0.37 ] 
0.12 - 0.13 10 0.39 - 0.40 3. 

649 
0.13 - 0.14 9 
0.14 - 0.15 6 
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TABLE F.2-7 


FREQUENCY DISTRIBUTION — OBJECT 64(S-band) 


o Interval No. of values o Interval No. of values 
(m*) in interval (m*) in interval 

0.00 - 0.01 49 0.27 - 0.28 6 
0.01 - 0.02 129 0.28 - 0.29 4 
0.02 - 0.03 147 0.29 - 0.30 3 
0.03 - 0.04 123 0.30 - 0.3] 6 
0.04 - 0.05 65 0.31 - 0.32 2 
0.05 - 0.06 52 0.32 - 0.33 2 
0.06 - 0.07 51 0.33 - 0.34 2 
0.07 - 0.08 59 0.34 - 0.35 3 
0.08 - 0.09 61 0.35 - 0.36 3 
0.09 - 0.10 73 0.36 - 0.37 3 
0.10 - 0.11 60 0.37 - 0.38 5 
0.11 - 0.12 28 0.38 - 0.39 4 
0.12 - 0.13 33 0.39 - 0.40 2 
0.13 - 0.14 19 0.41 - 0.42 ] 
0.14 - 0.15 15 0.42 - 0,43 5 
0.15 - 0.16 21 0.43 - 0,44 6 
0.16 - 0.17 14 0.44 - 0.45 6 
0.17 - 0.18 18 0.45 - 0.46 3 
0.18 - 0.19 18 0.46 - 0.47 4 
0.19 - 0.20 6 0.47 - 0.48 l 
0.20 - 0.21 11 0.48 - 0,49 7 
0.21 - 0.22 9 0.49 - 0.50 4 
0.22 - 0.23 7 0.50 - 0.5] 7 
0.23 - 0.24 8 0.51 - 0,52 3 
0.24 - 0.25 7 0.52 - 0.53 l 
0.25 - 0.26 6 1192 
0.26 - 0.27 10 
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o Interval 
(m*) 
0.00 - 0.01 
0.01 - 0.02 
0.02 - 0.03 
0.03 - 0.04 
0.04 - 0.05 
0.05 - 0.06 
0.06 - 0.07 
0.07 - 0.08 
0.08 - 0.09 
0.09 - 0.10 
0.10 - 0.11 
0.11 - 0.12 
0.12 - 0.13 
0.13 - 0.14 
0.14 - 0.15 
0.15 - 0,16 
0.16 - 0.17 
0.17 - 0.18 
0.18 - 0.19 
0.19 - 0.20 
0.20 - 0.21 
0.21 - 0.22 
0.22 - 0.23 
0.23 - 0.24 
0.24 - 0.25 
0.25 - 0.26 
0.26 - 0.27 
0.27 - 0.28 
0.28 - 0.29 
0.29 - 0.30 
0.30 - 0.31 
0.31 - 0.32 
0.32 - 0.33 
0.33 - 0.34 
0.34 - 0.35 
0.35 - 0.36 
0.36 - 0.37 
0.37 - 0.38 
0.38 - 0.39 
0.39 - 0.40 
0.40 - 0.41 
0.41 - 0.42 
0.42 - 0.43 
0.43 - 0.44 
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TABLE F.2-8 


FREQUENCY DISTRIBUTION — OBJECT 64(X-band) 


No. of values 
in interval 


56 
61 
54 
60 
52 
39 
46 
39 
33 
30 
4] 
26 
26 
21 
33 
27 
33 
35 
26 
25 
34 
27 
30 
34 
21 
16 
15 
10 
11 
13 
13 
12 
14 

6 
18 


o Interval 
(m”) 
0,44 - 0.45 
0.45 - 0.46 
0.46 - 0.47 
0.47 - 0.48 
0.48 - 0.49 
0.49 - 0.50 
0.50 - 0.51 
0.51 - 0.52 
0.52 - 0.53 
0.53 - 0.54 
0.54 - 0.55 
0.55 - 0.56 
0.57 - 0.58 
0.59 - 0.60 
0.60 - 0.61 
0.61 - 0.62 
0.63 - 0.64 
0.65 - 0.66 
0.66 - 0.67 
0.67 - 0.68 
0.69 - 0.70 
0.70 - 0.71 
0.74 - 0.75 
0.75 - 0.76 
0.80 - 0.81 
0.81 - 0.82 
0.83 - 0.84 
0.84 - 0.85 
0.85 - 0.86 
0.87 - 0.88 
0.88 - 0.89 
0.89 - 0.90 
0.94 - 0.95 
1.00 - 1.01 
1.04 - 1.05 
1.10 -1.11 
1.12 - 1.13 
1.13 - 1.14 
1.16 -1.17 
1.27 - 1.28 
1.28 - 1.29 
1.30 -1.31 
1.36 - 1.37 


No. of values 
in interval 
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TABLE F.2-9 


FREQUENCY DISTRIBUTION — OBJECT 65(S-band) 


o Interval No. of values og Interval No. of values 
(m4) in interval (m-) in interval 

0.00 - 0.01 78 0.20 - 0.21 ¢ 
0.01 - 0.02 44 0.21 - 0.22 5 
0.02 - 0.03 63 0.22 - 0.23 3 
0.03 - 0.04 74 0.23 - 0.24 6 
0.04 - 0.05 8] 0.24 - 0.25 7 
0.05 - 0.06 65 0.25 - 0.26 2 
0.06 - 0.07 75 0.26 - 0.27 5 
0.07 - 0.08 61 0.27 - 0.28 4 
0.08 - 0.09 44 0.28 - 0.29 6 
0.09 - 0.10 35 0.29 - 0.30 5 
0.10 - 0.11 3] 0.30 - 0.31 4 
0.11 - 0.12 38 0.31 - 0.32 5 
0.12 - 0.13 21 0.32 - 0.33 8 
0.13 - 0.14 31 0.33 - 0.34 5 
0.14 - 0.15 15 0.34 - 0.35 l 
0.15 - 0.16 9 0.35 - 0.36 2 
0.16 - 0.17 7 0.36 - 0.37 3 
0.17 - 0.18 14 0.37 - 0.38 S 
0.18 - 0.19 7 0.38 - 0.39 _— 
0.19 - 0,20 4 pee 


o Interval No. of values > Interval No. of values 
(m*) in interval (m*) in interval 
0.00 - 0.01 32 0.29 - 0.30 13 
0.01 - 0.02 49 0.30 - 0.31 9 
0.02 - 0.03 33 0.31 - 0.32 8 
0.03 - 0.04 33 0.32 - 0.33 5 
0.04 - 0.05 43 0.33 - 0.34 6 
0.05 - 0.06 51 0.34 - 0.35 10 
0.06 - 0.07 33 0.35 - 0.36 12 
0.07 - 0.08 38 0.36 - 0.37 7 
0.08 - 0.09 36 0.37 - 0.38 3 
0.09 - 0.10 31 0.38 - 0.39 9 
0.10 - 0.11 31 0.39 - 0.40 5 
0.11 - 0.12 33 0.40 - 0.41 l 
0.12 - 0,13 25 0.41 - 0.42 5 
0.13 - 0.14 27 0.42 - 0.43 9 
0.14 -0.15  . 24 0.43 - 0.44 6 
0.15 - 0.16 20 0.44 - 0,45 l 
0.16 - 0,17 17 0.45 - 0.46 Z 
0.17 - 0.18 23 0.46 - 0.47 4 
0.18 - 0.19 16 0.47 - 0.48 3 
0.19 - 0.20 19 0.48 - 0.49 7 
0.20 - 0.21 22 0.49 - 0.50 2 
0.21 - 0.22 16 0.50 - 0.51 l 
0.22 - 0.23 20 0.52 - 0.53 ] 
0.23 - 0.24 15 0.53 - 0.54 2 
0.24 - 0.25 23 0.55 - 0.56 2 
0.25 - 0.26 17 0.60 - 0.61 l 
0.26 - 0,27 12 0.62 - 0.63 l 
0.27 - 0.28 13 0.67 - 0.68 1 
0.28 - 0.29 10 898 
118 


5 [= CIRME TT 


UNIVERSITY OF MICHIGAN 
2428 -3-T 


TABLE F, 2-10 


FREQUENCY DISTRIBUTION — OBJECT 65(X-band) 
SS SSeS eather Pipe ress SSAA? 


S\|ECIRE IT 
UNIVERSITY OF MICHIGAN 
2428-3-T 


F.3 THE DEFINITIONS OF THE AUTOCORRELATION FUNCTION, 
R(T), AND THE POWER SPECTRUM, S(&), USED IN THE 
COMPUTATIONS 





In this section, a brief summary of pertinent facts concerning the 
usual R(t) and S(w), the definitions of R(T) and S(w) used in the computa- 
tions, the relationship between R(t) and S(w), and one aspect of the sig- 
nificance of the definitions are presented. 


When a real random function y(t) is known for all times t from 0 to 
co, the autocorrelation function RAT) and power spectrum, S(»), are 
defined as 





T 
R(t) = lim Ea | y(t) y(t+7)dt , (F. 3-1) 
Tp>o - 
0 
and 
2 
_ 2[A(w, T) | 
S_(v}) = lim ———— , F, 3-2 
oe ) T > T ( ) 
where 
T 
] -iut 
A(w, T) = | the at. F, 3-3 
ioe y(t) ( ) 
0 
Note that 
Si 20, = [Rr /<R,(0) , (F. 3-4) 
and 
Z 
S.(#) = s | R(t) coswrdr. (F. 3-5) 
0 


The last equation is known as the Wiener-Khinchine Theorem. 
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For a random function y(nTo), n= 0, 1, ..., N (that is, for a finite 
set of data), the autocorrelation function, R(t, T), and the power spec - 
trum, S(“), will be defined here by 


N-s 
] 
R(t,, T) = Nr : y(nT,) y(nTot+ts), OS TES NT, 
n=0 





= 0 , %7,<0 and T, PNT (F., 3-6) 
where 
Tg = STS, s = 0, 1, ..., N, and 
Z Z 
S(¥) = Alri, T) | : (F. 3-7) 
where 
N 
F(w, T) = ‘0 b, (nT ‘eo (F, 3-8) 
NT eg TNO | | 


Note that the summation in R(t,, T) is divided by N, the total number of 
values for y(nTo), rather than by N-s, the number of values of y(nT,>) x 
y(nTo+7g) in the summation. The reason for this is to preserve the 
relation 





N 
S(w) = ) R(sTo, T)coswsTy -5 R(O, T) >? , (F. 3-9) 


s=0 


which is the analogue of the Wiener-Khinchine relation. A proof of 
Equation F., 3-9 is the following: 


Using the definition of F(w, T) given in Equation F. 3-8, 
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Flu, T)|’ = 


2 N N 
O iw(n-m)T 
oD dD” ot) yemTgye(™T0 

n=0 m=0 


lusT 





0} y| (a -Ss)T ole 


n= s=n-N 





N N 
LZ PA y(nT ,) y|(n-8)T, |e'“*7 
s=0 ne=s | 


-] Nts 


? ey > y(n) y| (n-s)Tp Je#®T0 
s=-N n=0 





StS (nT,) +2 >» F wade s)T4| coswsT, 


s=l n=s 


From the definition of R(r,, T) 


N 
T y’ 

R(sTo, T) = (=e) y(nTo) y[(n-s)To] ; 
n=s 


where 
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Thus 


T 


N 

2 2 

S(%) = = Flo, 7)| = — R(o, T)+2 ) R(sT>, T)cos»sTy 
s=] 


from which Equation F. 3-9 immediately follows. 


One aspect of the significance of the definition of F(w,T) may be seen 
from the following considerations. R(7,, T) and S(w) were computed 
according to the above equations, 
using for y(nT,) the averaged 
values of cross-section de- 
f(t) 4 scribed in Section F. 2; that is, 
each o is the average over six 
: successive values of o with no 
overlap in neighboring t's. 
| Consider the Fourier analysis 
| of an arbitrary function f(t) 
| 
| 


muse =e OC 


sae aril which passes through a given 
set of points f,,, defined at the 
times a m=0, 1, ..., M, 
n=0, 1, ..., N, and f(t) = 0 
fort< le and t > ty, M+): (The 
Lay correspond to the o values, 
M+1 of these being used to ob- 


tain a value of ¢.) The Fourier 





t t > = 
oo Ol 36 6 


transform of f(t) is given by 


tn, M+1 N ™ ‘'n,m+! 
-luwt -iw 
F(w) = |—— f(t)e dt (a) S {(t)e dt , 
( =) | vém/ w=0 m= 
00 tom 
the latter expression being valid provided t, yqi) = ty4] 9 (this insures 


that the whole integration range be covered). 
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If 
tho mtl ~ ‘n,m = 4 anmdtog = 0, 
then 
tym = (Mtl) d+ma , 
and 
th, m+41 =n(Mt+l) A+ (mtl) aA , 
th, mtl A ~ivt! 
-1 -1 1 ! 
| tye it gt = oe i#tnm | dt! f(t'+tym)e 
thm 0 
If'wA<< 1, 
N M 
A -] es 
Flo) = ( ) Y _viw [n(M+1) Atma] ri, 
“T/n=0 m=0 
where 
n, mtl 
tnm 
If further wMA << 1, 
N 
F(o) = — s = “lunT'> 
‘nm n= 
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where 
= (M+1) A 


and 


M 
s ] eee 
'n > Ver Dy Tom - 


m=0 


The formal similarity between this expression for F() and f(w, T) as 
defined in Equation F.3-8 is apparent. If now f(t) is slowly varying over 
time intervals of length A (i.e., f(t) is band limited with “maximum 1/a), 
then fo ee t m+ 3 that is, fyy, may be replaced by fy, in which 
case f,, is merely the average of the Mtl values: fno. Ini. ---> fay 
Thus, “lie definition of F(w, T) given above and used in the computations 
gives the power spectrum, for w < 1/T 9 (To is the interval over which 
the data are averaged), of a function which is zero outside the observa- 
tion interval and which passes through all the points of the data before 
averaging. It is assumed that the function is band limited with maximum 
frequency << 1/A (A is the pulse repetition period of the radar used in 
obtaining the original data). In this case, A = 1/409.75 sec. and 

= 6A (M= 5). 


In actual practice, the following procedure was used: 


N-s 
R(t) = R(sT 9) = &. a, o(nT,) t ((nts) To) ’ 
n=0 


where values of s used were 


0(5)150(25)1250 for S-band , 


|] 
" 


s = 0(5)200(25)1250 for X-band . 


These are the values shown on the accompanying list. 
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Then S(w) for w< 2m was approximated as 





eT, x 
S(w) = ) ag R(sT,) cos(wST) , 
TT = 
where 
ag = 5/2 ; 
ae f 15<s< 145 for S-band 
a6." OF 15 <s< 195 for X-band ° 


150 for S-band 
= 200 for X-band ’ 


A) 
it 
foment 
mn 
bby 
oO 
ry 
—s\ 
nm M 
] WW 


oe. 175g <1250 for S-band 
a5. > OF} 225<s§<1250 for X-band 


This procedure was used to obtain S(w) for w/2m = 0(.1)1.0 cycles 
per second; for larger values of w, the method must be altered, inasmuch 
as taking s at intervals of 25 does not yield sufficient accuracy. In the 
interest of economy of time and money, it was decided to return to the 
basic definition of the power spectrum, namely 


S(w)= = P(e, T) i ; 


Numerically, this means 


IN/9 2 Tyg , 


o(5nT,) cos(5nuT o>) | + D>, F507) sin(5nwT 9) 
n=0 n=0 


25T 
mN 


e) 








S(¥) = 


This is the formula used to obtain S(w) for w/27 = 1.1, 1.2, 1.5, 2.0, 3.0, 
4.0, 5.5, 6.0, 6.5, 7.5 cycles per second. To get S(w) for higher values 
of w would require taking more than every fifth value of ¢ which was not 
considered warranted, 
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F,4 COMPUTATIONAL RESULTS FOR R(T) AND S(w) 


The results of the computations of the autocorrelation functions and 
spectra for the various objects are given in Tables F,4-1 — F.4-4 and 
Figures F,4-1—F,4-20, There are two curves for each object and fre- 
quency, one of which shows all of the computed values and is faired in at 
w/2n = 1,1 cycles per second, the other of which shows the region w<2z7 
in more detail. 


The spectra, as expected, all show a large peak at w=0. They then 
drop off rapidly and then oscillate more or less randomly. For only 
three of the ten cases did the ratio of the second largest peak to that at 
w = 0 exceed 1/10: for Object 59 at S-band this ratio is~0.2; for Object 
61 at S-band it is ~0. 3; and for Object 61 at X-band, there are two large 
subsidiary peaks with ratios to the zero-frequency peak of ~0.5 and 0.2. 
The observation times for Objects 59 and 61 were considerably less than 
for the others. (The sample lengths were approximately 7, 3, and 3 
seconds for Object 59 at S-band, Object 61 at S-band, and Object 61 at 
X-band, respectively. For the other objects, the sample length was >10 
seconds. ) 


Before interpreting the peaks in terms of physical characteristics 

of the missile, one should take into account the fact that the zero- 
frequency signal from the missile will, for short observation times, con- 
tribute appreciable subsidiary peaks to the spectrum (in addition to the 
one at w= 0). When this is done, it is found that the results are not in- 
consistent with the assumption that the spectra are due merely to a zero- 
frequency signal with superimposed noise characteristics inherent to the 
small sample. If the missile generates any other frequencies, these are 


not likely to be detected from the power spectra. 


F.5 FURTHER ANALYSIS OF R(T) AND S(w) 


The finite nature of the sample involved in these computations of the 
autocorrelation function and power spectrum for the radar cross-section 
data obtained in the drop tests necessitates a more detailed study, which 
is presented in this section. 
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15 
20 
25 
30 
35 
40 
45 
50 
55 
60 

65 
70 
75 
80 
85 
90 
95 

100 

105 

110 

115 

120 

125 

130 

135 

140 

145 
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Object 59 
-, 0008453 


0008353 
0008000 
0007538 
0006946 
0006308 
0005674 
0005087 
0004469 
0003924 
0003453 
0003045 
0002761 
0002554 
0002460 
0002395 
0002451 
0002522 
0002656 
0002725 
0002775 
0002768 
0002732 
0002676 
0002540 
0002395 
0002230 
0002076 
0001864 
0001721 
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TABLE F.4-1 


ag 


R(n T,) —5-BAND — SWEPT BACK FINS 
Object 61 
- 0029390 


0023537 
0017303 
0011459 
0007092 
0005899 
0005798 
0006064 
0007060 
0007693 
0008243 
0009202 
0010578 
0012257 
0014454 
0013674 
0011573 
0009491 
0006780 
0005560 
0005202 
0004491 
0003977 
0003463 
0004005 
0005413 
0006821 
0008294 
0009101 
0009339 


Object 62 
. 0023321 


0018821 
0013379 
0009538 
0007803 
0007903 
0008497 
0009763 
0010370 
0010078 
0009771 
0010260 
0010491 
0010660 
0010257 
0009476 
0008434 
0007511 
0007332 
0007779 
0008361 
0008304 
0007481 
0007005 
0007566 
0008530 
0009498 
0010403 
0010378 
0009784 


150 
175 
200 
225 
250 
275 
300 
325 
350 
375 
400 
425 
450 
475 
500 
525 
550 
575 
600 
625 
650 
675 
700 
725 
750 
775 
800 
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TABLE F.4-1 (Cont. ) 


0001308 
0001248 
0001098 
0001533 
0001949 
0001547 
0000730 
0000348 
0000442 
0000435 
0000129 
0000000 


Sle CIRIE 


128 


R(n T,) —S-BAND — SWEPT BACK FINS 
Object 59 
.0001623 


Object 61 
0007954 


0005069 
0002642 
0006688 
0000477 
0000812 
0003202 
0000193 
0000615 


MICHIGAN 


Object 62 
.0008917 


0006795 
0007798 
0006138 
0006063 
0007759 
0007138 
0005265 
0007893 
0005356 
0003522 
0006118 
0003889 
0003649 
0003116 
0002660 
0002799 
0002389 
0003582 
0003334 
0002809 
0002497 
0001704 
0002867 
0000575 
0002067 
0000082 
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TABLE F.4-2 


R(n To) — X-BAND — SWEPT BACK FINS 


n Object 59 Object 61 Object 62 

0 0141270 . 0046145 .0050741 

5 0071525 0041448 0025862 
10 0070421 0034787 0020497 
15 0059138 0028063 0017938 
20 0053559 0021480 0021512 
25 0044026 0015557 0025621 
30 0054613 0010679 0021173 
35 0043404 00060991 0018398 
40 0040270 0004715 0018803 
45 0045957 0003385 0023149 
50 0041832 0002416 0022695 
55 0033812 0001878 0016563 
60 0034853 0001489 0018283 
65 0033959 0001520 0018065 
70 0032140 0002208 0019548 
75 0034870 0003448 0022300 
80 0036184 0004371 0017712 
85 0036834 0005240 0021593 
90 0044616 0005950 0017755 
95 0055883 0007081 0017029 
100 0056242 0008855 0016308 
105 0056562 0011552 0018463 
110 0049853 0014348 0017689 
115 0041834 0016783 0017517 
120 0039607 0018905 0015197 
125 0036294 0020430 0015692 
130 0037868 0022077 0019122 
135 0037045 0022308 0018591 
140 0036235 0020588 0017432 
145 0033121 0017566 0017085 
150 0029456 0014172 0014481 
155 0026168 0009584 0014074 
160 0024110 0005271 0012283 
165 0024711 0002434 0010523 
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170 
175 
180 
185 
190 
195 


200 
225 
250 
275 
300 
325 
350 
375 
400 
425 
450 
475 
500 
525 
550 
575 
600 
625 
650 
675 
700 
725 
750 
775 
800 
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TABLE F.4-2 (Cont. ) 


Object 59 
. 0026162 


0026965 
0027853 
0029782 
0032698 
0035881 


. 0039490 


0034259 
0022680 
0021086 
0022799 
0023015 
0013603 
0012445 
0009518 
0007890 
0002006 


R(n To) ~ X-BAND — SWEPT BACK FINS 
Object 61 
. 0001081 


0000774 
0000973 
0001878 
0003357 
0003873 


. 0003661 


0003475 
0002882 
0001199 
0000566 
0001348 
0001063 


Object 62 
0010752 


0012779 
0013834 
0011650 
0013454 
0010824 


.0011133 


0011588 
0009446 
0007613 
0006469 
0008158 
0006938 
0005678 
0004964 
0003596 
0004875 
0004372 
0004175 
0004257 
0004977 
0003584 
0005147 
0003615 
0002293 
0002339 
0001387 
0000443 
0000183 
0000331 
0000000 
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10 
15 
20 
25 
30 
35 
40 
45 
50 
55 
60 
65 
70 
75 
80 
85 
90 
95 
100 
105 
110 
115 
120 
125 
130 
135 
140 
145 
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TABLE F.4-3 


R(n Ty) — S-BAND — SQUARE FINS 


Object 64 


0213705 


S55 


0200862 
0178543 
0156342 
0138965 
0125703 
0115042 
0106352 
0098932 
0092760 
0088143 
0085914 
0087777 
0091778 
0095362 
0096523 
0095611 
0093766 
0091029 
0087031 
0082757 
0078319 
0073477 
0068420 
0064605 
0063129 
0063519 
0064878 
0067283 
0071034 
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Object 65 
. 0135879 


0102389 
0066842 
0061679 
0063365 
0068854 
0075865 
0079173 
0077494 
0064617 
0057250 
0060818 
0065239 
0066173 
0067007 
0067176 
0061667 
0054871 
0055849 
0063276 
0073871 
0073213 
0058965 
0049297 
0047688 
0051110 
0058077 
0059460 
0058439 
0057026 
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TABLE F.4-3 (Cont. ) 


R(n T,) — S-BAND — SQUARE FINS 


n Object 64 Object 65 
150 . 0075680 . 0052668 
175 0089139 0053268 
200 0077239 0059709 
225 0060214 0045687 
250 0063744 0037299 
275 0066082 0042200 
300 0045265 0044469 
36) 0061719 0036720 
350 0070511 0036091 
375 0056870 0037898 
400 0045454 0040035 
425 0045930 0028989 
450 0041060 0029714 
475 0040726 0032119 
500 0048994 0039226 
525 0048728 0029245 
550 0043212 0025457 
575 0062451 0028174 
600 0067761 0038098 
625 0050572 0022803 
650 0062977 0026753 
675 0072518 0038686 
700 0056768 0028555 
725 0039070 0018952 
750 0035771 0017655 
775 0032301 0015259 
800 0027729 0009849 
825 0026783 0007866 
850 0051205 0007180 
875 0035928 0008692 
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TABLE F.4-3 (Cont. ) 


R(n To) —-S-BAND — SQUARE FINS 


n Object 64 Object 65 
900 0019334 . 0004677 
925 0036836 0001206 
950 0027423 
975 0013044 

1000 0020053 
1025 0023112 
1050 0016797 
1075 0007424 
1100 0005445 
1125 0009259 
1150 0002440 
1175 0002134 
1200 0006488 
1225 0006437 
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TABLE F., 4-4 


R(n Tg) — X-BAND — SQUARE FINS 


n Object 64 Object 65 

0 .0731492 . 0401553 

5 0496405 0269065 
10 0415901 0257876 
15 0418286 0257512 
20 0359315 0245474 
25 0301876 0238165 
30 0310012 0244880 
35 0311995 0242412 
40 0350713 0226309 
45 0377331 0255553 
50 0422809 0256857 
55 0432972 0253318 
60 0370850 0232201 
65 0359330 0252210 
70 0362422 0235390 
75 0348520 0251941 
80 0299206 0224900 
85 0289381 0217394 
90 0330233 0238443 
95 0332032 0240687 
100 0329484 0243688 
105 0400315 0241638 
110 0364368 0232737 
115 0347348 0216312 
120 0358902 0237124 
125 0335385 0218293 
130 0290651 0205318 
135 0272524 0213441 
140 0280468 0226513 
145 0301092 0235193 
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TABLE F.4-4 (Cont. ) 


R(n T,) — X-BAND — SQUARE FINS 


n Object 64 Object 65 
150 . 0291169 .0237621 
155 0333597 0222940 
160 0355673 0209540 
165 0326777 0225274 
170 0333055 0221028 
175 0351747 0212266 
180 0319484 0205320 
185 0296904 0199072 
190 0284441 0206678 
195 0282039 0225460 
200 0279145 . 0220419 
225 0300524 0190359 
250 0251444 0181058 
275 0272144 0174472 
300 0260056 0159126 
325 0256764 0159699 
350 0250086 0152185 
375 0246935 0145686 
400 0261566 0139706 
425 0237105 0112744 
450 0274815 0114447 
475 0204867 0109566 
500 0229171 0098070 
525 0226876 0098986 
550 0204659 0084603 
575 0185845 0081030 
600 0181684 0063576 
625 0173128 0052469 
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650 
675 
700 
725 
750 
775 
800 
825 
850 
875 
900 
925 
950 
975 
1000 
1025 
1050 
1075 
1100 
1125 
1150 
1175 
1200 
1225 
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TABLE F.4-4 (Cont. ) 


R(n Tg) — X-BAND — SQUARE FINS 
Object 64 


0151293 


5 |= 


0148216 
0137398 
0133801 
0088473 
0089325 
0077180 
0069846 
0080214 
0073055 
0052689 
0053285 
0040715 
0039828 
0035806 
0026802 
0031196 
0031564 
0024238 
0020768 
0019607 
0014503 
0011415 
0003088 
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Object 65 
. 0049851 


0041182 
0037598 
0032464 
0021562 
0017225 
0013484 
0008427 
0007767 
0008462 
0003787 
0002023 


S (w) x 104 (meters 4 / cps) 
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F,5.1 Definition of the Autocorrelation Function and Power Spectrum 
for a Continuous Sample on a Finite Interval 


In defining the spectral density (power spectrum) and autocorrelation 
function when the data is given for only a finite interval (2T), these quan- 
tities must approach, continuously, the corresponding definitions given 
in Equations F. 3-1 and F. 3-2 as T becomes larger. Definitions for the 
continuum corresponding to those of Section F. 3 are 





5,(4,T) = a ne (F.5-1) 
] 
Ry(7,T) = +: | y(t, T) y(te7, T) dt, (F. 5-2) 
where 
T 
l -iut 
F,(¥,T) = the dt , 
-T 
y(t), Itl< T 
y(t, T) = : (F. 5-3) 
0, It] >T 


These functions have the properties S1(,T) 20, Ry(7, T) = Ry(-7,T), 
R}(0,T)2 Ry(7,T), and they Satisfy the Wiener-Khinchine relation: 


Z 
91 (4, T) = ve 


oO 
| Ri(7,T) coswtdr , (F. 5-4) 
0 
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The proof of Equation F.5-4: By definition, 


3 
8 


l oe ae 
91(4, T) = oT dt dt’ y(t, T) y(t', T)e TNE 


i 
8 
I 
re) 


i r : 
LWT 
-3T | dt y(t, T) dt y(t+7, T)e 
=-0O -~& 
aad ° 
=-— | dre” Ry(r,T) . (F.5-5) 


From Equation F’. 5-2, it follows that R)(-7,T) = Rj(7,T), so that Equa- 
tion F.5-5 implies Equation F.5-4. Inverting, we have 


oO 


; 7 
Ry(7,T) = > | dwe”” $i(w,T) . (F. 5-6) 
—0o 
Thus, 
] OO 

Run TI<> ( dw Sy(u,7) , 
or 

Ry (7, T)| <R,(0,T) . (F. 5-7) 
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Since the limiting value, R(r), is the time average of the product 
y(t) y(t+r), another reasonable definition, R2(rt, T), would be 


l r 
Rats T) = ae | y(t, T) y(tt7, T) dt ; (F. 5-8) 


that is, the integral 


ore) T-T 
| y(t, T) y(tt+7, T) dt = | y(t) y(t+7) dt, 7 20 
=O -T 


is divided by the interval used in averaging y(t) y(tt7). Note that 
\R2(7,T)| can be greater than R,(0,T). For example, if y(t) is as shown 
in the figure, R2(0,T) = tp/T, whereas R>(2T-to, T) = 1 > R2(0, T) for 
0<t,< T. 





One could also define a spectral density as the transform of R(t, T): 


oO 


2 
S2(#, T) = = | R2(7,T) coswtdr, (F, 5-9) 
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That 52(w, T) is not non-negative may be seen by the example y(t) = 1. 
Then R2(7,T) = 1, and S9(w, T) = (2/7) sin 2wT. Note that 


2T 


ae ae Cs 
51(4, T) = | (: =) R2(7t,T) coSwrdr , 
0 


1,€., 5j(w,) is the first order Césaro sum of S?(w, 0), so that when the 
latter exists it is equal to Sj(u, 0). 


Before comparing the relative merits of these various definitions, 

the question as to what is desirable in the approximate quantities must 
be discussed, Clearly the S(w, T) and R(t, T) which minimize IS(w, T)-S(w) 
and | R(7, T)-R(7)| , for all w and +, respectively, would be the most de- 
sirable. On the other hand, the usual types of criteria, wherein a single 
parameter such as 

[ 

0 


ta { fst.m-s00)] a 
0 





Ly S(w, T)-S(w) dw 








or 


is minimized would not necessarily yield the most desirable approximate 
function, For example, if y(t) consists of a random (noise) component 
and a signal made up of a finite number of sine waves, and the purpose 

is to detect this signal and determine the frequencies, one would choose 
the S(w, T) which gave the "sharpest'' peaks at (or near) those frequencies, 
provided it did not give sharp peaks at other frequencies. It can be 
Shown that, in general, neither L) nor L> are good comparison param- 
eters for an analysis of this kind. 
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Criteria involving S(v) and R(r), the true spectrum and correlation 
function, cannot be applied directly to problems such as the data analysis 
for the Bendix Aviation Corporation, since the limitation to a finite 
sample excludes the possibility of knowing S(w) and R(t). However, it is 
conceivable that some process of guessing the statistics and finding opti- 
mum definitions for R(t, T) and S(w, T) could lead one to better definitions 
(Ref. 29)". This approach, which seems worthy of investigation, will 
not be considered further here. 


‘In Reference 29, the finite time average, 
T 


| Z(t) dt, 


where Z(t) is some functional of y(t) (both y(t) and Z(t) are assumed to be 
stationary) is generalized to 
ah 


| g(t) Z(t) dt =M(T); the quantity <r) - <r) >| ‘un? = 
0 


Zs 
T 


is then minimized ({ }denotes an ensemble average). This leads to an 
integral equation for the "optimum" g(t), having the autocorrelation func- 
tion <Z(t) Z(tt7)> as its kernel. One can apply this to the autocorrela- 
tion function <y(t) y(t+7)> by taking Z(t) = y(t) y(tt7) and, in addition, 
modifying the limits of integration to 

T-7T 

| 
However, the process of deriving the spectrum from the optimum auto- 
correlation function through the Wiener-Khinchine relation does not, in 
general, give an optimum definition for the spectrum. Further, the 
analysis of Reference 29 is of no help in trying to minimize L directly 
for the spectrum, i.e., taking 

T 
Z(t) = { K(t) y(t)e” at. 


= E 
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Other properties which might be considered desirable are those 
described by Equation F. 3-4. These will not be assumed here (if they 
were, then the pair R2(7, T) and S2(w, T) would immediately be ruled out). 
Instead, y(t) = Acos wot + N(t) will be investigated, where the function 
N(t) is assumed to have some standard noise properties. S(w,T) will be 
required to consist of a continuous (noise) component plus a peak’ located 
near w= wo(the purpose of the peak is to allow one to say whether or 
not the signal, AcoSwot, is present and, further, to determine the value 
of wo). This will be done for the two functions S)(w, T) and S2(w, T); it 
will be seen that the nature of the signal peak, and in fact the whole func- 
tion, is essentially different for these two definitions. Some remarks 
about the comparison of Rj(7, T) and R2(7, T) will also be made. 


For y(t) of the form 
y(t) = o(t) + n(t) (F. 5-10) 


(o(t) is considered to be signal, n(t) to be noise), 


Rj(7, T) = = | U(t, T) U(ttr, T) ott + n(t) Letts) + n(te7)| ; 


where 


‘It is tempting to consider the peaked component an approximation to 
the Dirac $-function. Care should be taken in this respect, however, 
since there the manner in which the approximation approaches the 
é6-function is important. For example, the function 8(u, T) = J 2/n Tcosw“T% 
satisfies the integral property 


J bw-0" ful) Got Sku 


in the limit T>~for most well-behaved functions f(x), but this 4(w, T) 
would not be satisfactory for our purposes. 
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Thus, 





T- 
- Ro n aun : 
R(t, T) = R (7, T) + R (7,7) + aT [| o(t) n(t+7) + o(t+7) n(t) Jat, 
2T>7Td0, (F.5-11) 
where 
T- 
R Mr, T= = [ B(t) B(t+r7) dt. B=a,n. 


Taking the ensemble average, 


T- 
R (7,7) = R ee T) + dp" (Zz; T)> + Sar sl an ote + o(t+7) at, 


(S77 . (F, 5-12) 


where n(t) is assumed to be a stationary random process; this assump- 
tion also leads to 


<Ri (7, T)> = (1-=) R" (7), O<7T<2T , (F. 5-13) 
where 


R"(r) = n(t) n(ttr)> (F. 5-14) 


Noting that 





R,(7,T) = 


R(t, T) . (F. 5-15) 
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and 
raah 
B _2 B 
5, (o, T) aoe | R, (7, T) coswtdrT , 
0 
then 
oT T-T 
<S, (w, T)> = S ile, T) + <si (w, T) +2 | ar cosur | oft) + o(t+7) | dat 
0 -T 
Zy: 
=§ ila, T) + Csiu, T)> are | dr sinur|o(T-r) to(r-T) ; 
0 
(F. 5-16) 
and 
aT. T-+ 
Z WT 
‘S, (w, T)> = S stu, T) + <8) (w, T)> + 20 a dr | dt fo(t) + o(t+>) | : 
0 -T 
(F, 5-17) 
Clearly 
or 
(Su T)> = e | 0-) R" (rT) coSwtdT (F. 5-18) 
hea T 2T : 
0 
and 
2T yaa 
(S3(w, T)> -< | R"(r) COSwTdT = (Si (u, T)»> — | Rn coSwtdr . 
0) 0 


(F. 5-19) 
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Assuming <1>= 0, the autocorrelation functions and spectra are simply 
sums of the contributions due to signal alone and noise alone. Further- 
more, from Equation F. 5-19 it follows that for T >>the correlation time 
of the noise, 


(Siw, T)) = CSh(w, T)> . 


Consider the part of the correlation function and spectrum due to signal 
alone. For 








a(t) = AcoSuwot , (F. 5-20) 
simple integration leads to 
ss onli 1 Sinw(2T-T) nee 
R (7, T) =A AT coswgr Pot a , OST <Z27 . (F. 5-21) 
v nae 1 sinw(2T-T) oe 
Ri(7, T) =A 3 cos WoT + Tor a 4 USAT (F’, 5-22) 


ae a2 sin"(o-v)T sin@(utu,)T _ 2 Sinle-wolT sinlotug)T 
tine 2nT (w- wo)? (aeaey" (w-Wo) (wt) 


(F. 5-23) 
Calculating S5(w, T) from the definition, 
2}/sin2(w-w.)T  sin2(w+w,)T 2 
A 
Solu, T) = —— A 


2m +7 ltv,T), (F. 5-24) 
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where 
raat 
COS WT Sin w(2T-T) 
I(w, T) = dt ITo7 | 
0 


It is easy to show that 


I(w, T) = cos ieT {si [2(wt+wo) T] ~ si(zo-uo) + — {v te-)a}} P 








(F. 5-25) 
where 
x 
, sint 
51(x) | ; dt , 
0 
x 
l-cos 
W(x) = | ; : dt = Inyx-Ci(x) , 
0 
ie ©) 
Ci(x) = —— dt , and 
x 
Y = 1,781072, Euler's constant. 
Consider first the spectra for the case Wy = 0. Then 
2 bcc 
S"(u, T) = sa ed (F, 5-26) 
and (uw, = 0) 


sw, 7 _ 2A* sin2wT 
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(since I(w, T)/w = (sin2wT)/w when W) = 0). The interesting feature here 
is that, while for T+>- both functions satisfy, to within a normalization 
constant, the property 


| (0) f(w) dw = £(0) 


(for most well-behaved f(w)), they differ in a number of other essential 
properties. Some of these are advantageous (from the viewpoint of 

Signal detection and resolution) for one, some for the other. For ex- 
ample, the fact that S3(w, T) is twice as high as Si(w, T) at w= 0, and 

half as wide (i.e., the first zero of So(w, T) occurs at half that of S1(w, T)), 
works in favor of S5(w, T). On the other hand, the property So (4, T)TS OQ, 
“= constant # 0, which is not exhibited by S5(w,T), tends to favor Su, T). 
Another property which favors Si(w, T) is that the ratio of the height of 

the peak at w = 0 to that of the neighboring peak (which occurs at 

W = W} = 37/2T) is almost three times the corresponding ratio for S3(w, T): 





o cv, , 
$)(0, T)/S (4). 7) 9n*/4 On | 
$$ se =—~-=2,8 , (F. 5-27) 
s(0,T)/s%w!,T) 7%} 
(0, T)/S5(o!, T) 
Consider w) >0. Then, 
5° A2T sin*(w-,)T 
( T) sae ta (Rear oma A \(w, T) ; (F. 5-28) 
/ Aé@T |Sin2(w-w5)T 
5,(4, T) wALT (tee : (F, 5-29) 
where the functions 
Sin(w-w5)T sin(wt+w )T sin*(w+w )T 
A,(w,T) = 2 = & ———— , (F. 5-30) 
(w - Wo)T (w+wWo)T (w+) Te 
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Sind(wt+w,)T I(w, T) 


(F. 5-31) 
tend to distort the peaks; these peaks are represented by the terms 


(A*/2nT) [sin2(o-u9)T | / (w-wWo)° 


and 


(A“/2n) [sin2(o-uo)T | / (we) 


in Si lw, T) and $5(u, T), respectively, which are simply 1/4 the functions 
for wy = 0 (Eq. F.5-26), but centered around w>. Clearly there is not 
much distortion in Stl, T) for w™) T>>z; for example, Aj(w,T) = O(w 7! To )) 
at w= 4. To estimate the distortion in S5(u, T), examine the asymptotic 
expressions (Ref. 30) 


T  CcOSx sin x 
Si(x) © > - — ’ ¥(x) = Iny x - P 


valid for x > 2m. Considering again the value at w= u4,, 


Ssin4wo5T 


—__—_——_—_— ~ 5-32 
4uoT : (F 


Aol, T) = [in 4> WoT) sin2 wT + cos 2u9T | + 


which goes to zero as uw) T>~; the logarithmic term, however, tends to 
cause more distortion in the case of S5(w, T). 


Thus the two functions, 5](w, T) and S2(w, T), which differ in several 
essential respects, both show peaks at the signal frequency, plus other 
maxima which are less than 1/7 the height of the signal peak. For the 
purpose of the Bendix contract, wherein the only a priori information 
was that there was a large zero-frequency component, with the object 
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being to discover any other components if they existed, the choice ve- 
tween the various definitions was a matter of indifference. 


Before concluding this section, one further remark concerning the 
two definitions of the autocorrelation functions remains to be made. 
R,(7, T) is identically zero for +22T and is near zero for 7=2T, where- 
as R2(t, T) will not generally approach zero ast>2T (7 <2T). There- 
fore,R2(t, T) could be considered, in general, a better approximation to 
the true R(r). However, as t>2T, the length of sample used in com- 
puting either R(t, T) goes to zero, so that the results for 7r=2T will be 
statistically unreliable; hence, essentially nothing is lost in this respect 
by using R}(7, T). 


F.5,.2 The Autocorrelation Function and Power Spectrum for a 
Finite Sample 

The drop-test data, being the result of measurement with a pulsed 
radar (prr =400 sec”), yielded a finite number of values of the cross- 
section spaced at intervals A=1/400 sec, The problem of defining 
the autocorrelation function and power spectrum in terms of such data 
is essentially the same as in the case of a continuous sample of finite 
length with the exception that, in the discrete case, one must make an 
assumption as to band limitation. This corresponds to the fact that 
cos(w ) t) and cos(w2 t) will look precisely the same if sampled at times 
th =nA(n=0, 1, ..., N) if #2 = & + 2s7/A, S any integer. The limita- 
tion imposed by the discrete observations is not restricted, a priori, 
to low frequencies. It is necessary only that two frequencies differing 
by 2s7/A be indistinguishable (no matter how large N may be). However, 
physical considerations lead to the assumption that the function being 
sampled is essentially band-limited by w<27/A =400 cps (that is, the 
amplitudes of the higher frequency components are much less than the 
amplitudes of those within this band). The actual computation was not 
based on the original data, but on a reduced set obtained from the origi- 
nal by averaging, with no overlap, successive sets of six points. Hence, 
the data used was at intervals of Tp = 64. One may ask whether it is 
necessary to assume a band limitation of w<27/Tp = 66 cps. It was 
shown in Section F.3 that, in terms of the definitions used, this was not 
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the case — one need assume only the band limitation imposed by the 
spacing of the original data, provided that the spectrum is computed for 
w <<2n/To. 


Analogously to the case of a continuous sample, there are, as alter- 
natives, the definitions of Section F. 3: 


: N-s 
] 
R'(t5,.T) = Way y ; y(nT,) y(nT9 + Ts) 


and 


N 

2 l 

S'(w) = “2 S R'(sTo, T) coswsTo os R'(0, T) 
s=0 


For the case y(t) = 1, the summations are easily carried out, yielding 


sin (N+l)wTo 


T 2 
sin@ v0. 
S'(w) =—— <———— + cosNoT, 
‘ WloO 


Hence, the remarks of Section F.5.1 concerning the comparison of the 
various definitions apply here as well. 
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APPENDIX G 


THE CROSS-SECTION OF THE TANKS 


G.1 INTRODUCTION 


The cross-section of the terminal propulsion unit (tanks) at 200 Mc 
was determined by essentially three different methods of approach: 


1, Using previously obtained theoretical and experimental 
data on the 7-OC booster (Ref. 1), 


2. Using physical optics methods on the breakdown of the 
tanks shown in Figure G.1-1 (a drawing of the tanks 
appears in Fig. 2-4"), and 


3. Using physical optics methods on a modification of the 
breakdown used in Method 2. 


In the consideration of the problem by Methods 2 and 3, the cross- 
section was determined for:the tanks with the warhead attached, with 
the warhead removed, and with both the warhead and hat section removed. 


G.2 METHOD 1 —SCALED 7-OC BOOSTER DATA 


This section contains an estimate of the radar cross-section of the 
terminal propulsion unit (tanks) at 200 Mc based upon the results obtained 
previously on the 7-OC booster. Experimental and theoretical data on 
the cross-section of the 7-OC booster appear in Reference 1. If the 
7-OC booster is scaled down by a factor of 5/6 (Fig. G.2-1), the result- 
ing model closely resembles that of the tanks. The dimension of the 
5 /6-scale 7-OC booster and the tanks are compared in Table G. 2-1. 


'The drawing of the tanks shown in Figure 2-4 shows the warhead in 
the Shape of a cone rather than as one of the Rudolph parameterizations. 
Thus, the Method 2 approach, in computing the cross-section with the 
warhead attached, assumes the warhead to be in the shape of the cone 
Shown on Figure 2-4, 
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Hat Section 





ey 


Nose Cone 


FRONT SECTION 


CENTER SECTION 









(Ogive: P=100") 
50" REAR SECTION 


a — cylinder: 103” long; dia. = 10” (sphere cap) 
b — cylinder: 50” long; dia. = 18” 

c — sphere: dia. = 15” 

d — sphere: dia. = 20” 


FIG. G.l-1 GEOMETRY OF “TANKS” USED IN COMPUTATIONS 
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(all dimensions are in inches) 


FIG. G.2-1 A 5/6 -SCALE MODEL OF THE 7-OC BOOSTER 


TABLE G. 2-1 
COMPARISON OF DIMENSIONS BETWEEN THE TANKS AND THE 


5 /6-SCALE MODEL OF THE 7-OC BOOSTER 
7-OC Booster 


Maximum diameter 
Front diameter 

Rear diameter 
Length-truncated cone 
Length-cylinder section 


Length-rear sloped section 


Length-tail cylinder 


Nose angle (truncated cone) 
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Application of modeling theory to the data for the 7-OC booster 
given in Reference 1 yields cross-section data on the '5/6-scale 7-OC 
booster" (and thus on the tanks) at 90, 270, 720, and 1200 Mc. Figure 
G.2-2 shows the theoretical cross-section (average oc) of the "5/6-scale 
7-OC booster" at 270, 720, and 1200 Mc,as scaled from the data of 
Reference 1. 


The theoretical curves given in Figure G. 2-2 for the cross-section, 
combined with the experimental data shown in Figures G.2-3 and G.2-4', 
yield a reasonably good estimate of the cross-section pattern to be 
expected for the tanks at 200 Mc in the interval 0< ¥<90°. 


G.3 METHOD 2 — PHYSICAL OPTICS CROSS-SECTION OF THE TANKS 





G.3.1 Cross-Sections 


The physical optics method was applied to the determination of the 
cross-sections of the tanks using the formulas given in Section G.3.2 
which were evolved from the breakdown of the tanks by the methods of 
Appendix A of Reference 5. The cross-sections were determined at 
200 Mc for the tanks with: 


1. The nose-cone attached; 


2. The nose-cone removed and the resulting ho.e assumed 
to be: 


a. Filled with an absorbing material so that no 
return is received from the interior of the hole, 


b. Lined on the inside with an inverted hemisphere, 
and | 


c. Empty up to the end of the hat section, so that 
from the nose-on aspect the return would be 
like that from a flat plate; and 


3. Both the nose-cone and hat section removed. 


‘Figures G.2-2 —G.2-4 were taken directly from Reference 1, with 
the appropriate scale factors introduced. 
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FIG. G.2-2. THEORETICAL CROSS-SECTIONS OF THE TANKS 


(Based upon a 5/6-scale model of the 7-OC Booster) 
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FIG. G.2-3 CROSS-SECTION OF THE TANKS (NOSE-CONE REMOVED) AT 90 Mc 


(Based upon a 5/6-Scale Model of the 7-OC Booster) 
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(Based upon a 5/6-scale model of the 7-OC Booster) 
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The results obtained are shown in Figure G.3-1. 


G.3.2 Computational Procedures For the Tanks 


In determining the cross-section, using the breakdown shown in 
Figure G.1-1, the primary contributions to the cross-section were 
assumed to come from: | 


1. The tip of the nose-cone (cj), 
2. The junction of the nose-cone and the hat section (o, and v,), 


3. The junction of the extended hat section and the central 
cylindrical portion («7 and os), 


4. The junction of the cylindrical portion and the rear 
conical surface (og and og), and 


5. The engine breakdown back of the rear cone (c7). 


At certain aspects (namely those which are normal to the conical or 
cylindrical surfaces), other formulations were involved. 


The formulas involved are listed below: 


2 2 2 2 
ee R,, tan (10°) +7R, [tana - tan (10°) , for y= 0° 
5 
O O 
o = p v. », tor 0 <w<l10 
i=] 
6 
o =o, + 3 o. , for 10°<¥< a 
l 1 
=3 
7 
c=o_, +t ae , for a<v<909 -a 
] co4 1 
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8 
— ) c, , for 170°<y< 180° 
i=6 


2 
g = n(r4)° +(R2)* tan a+o7 , for ¥= 180° 
where 


A=1.5m 
a = tan7! (24/29) 
= 0.61 m 
R,, = 1.52 m 
r=0Q.23m 
r= 0.63 m 


and 


3 

7 * tan*a 1+ cos(2a) 

l lot cos(2¥) +cos(2a) 
AR, tan“(10° - y) 

2 87 siny 


AR 


2 
_ oe oO” O O 
v. = irony (tan (a - ») - tan (10 | » ford <¥<10Q 


_ AR] 2 O 
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AR, ; Z 
v, = Tea tent19 + v) 7 tans) } 


AR> 2 5 
°° semi v} , for 10 <wca 


AR? a 4 
= _——__¢ |tany| - |tan(p - a)| , for a< p< 180° - 
8n | siny| a a 
(yp 7 a oy + a) 


2 
o => (51) + (.37) fs fora <p<90° 


N| A 


2 2 2 
(51) + (.37) + (.25) | , for 90° << 180° 


Z 
= (37° + (.25) , for » = 180° 


Ro 
T = maunin” (a + a} 


The formulas above were applied for the tanks with the nose-cone 
attached. In considering the cases of the tanks with the nose-cone 
removed, the expEcoe=ons for oc, i= 1, 3, and 4 had to be changed and, 
at the aspect ¥ = 90° - a, the normal incidence term deleted. 


The new forms of v, and c, are 


AR) 
! ee ; 
3 St Sin 


! 


tan“ (10° - »), 0°<v< 10° 





AR, 
1 
T= 
4 8m siny 





tan“ (10° + vw), wF g0° ; 
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and v) is replaced by: 


1. The wire loop formula for the case cited in (2a) of Section 
G.3.1, 


2, Vhe wire loop and sphere formulas for the case cited in 
(2b), and 


3. The wire loop and circular flat plate formulas for the 
case cited in (2c). 


In determining the cross-section for the tanks with both warhead 
and hat section removed, the approach outlined above for the case cited 
in (2c) was employed, with the value of R; changed to 0.88. 





G.4 METHOD 3 — CROSS-SECTION OF MODIFIED TANKS 


In performing the analysis and computations for the Method 2 
approach, it was observed that the base of the nose-cone shown on the 
drawing furnished to us of the terminal propulsion unit was larger in 
diameter than the maximum diameter of the Rudolph warhead. Thus, the 
tanks were modified slightly (the hat section was lengthened), and the 
computations were repeated. The warhead used in these computations 
is Configuration IV. 


Configuration IV was chosen because it most closely resembled the 
nose-cone shown on the "tanks" drawing. One further modification was 
made in the breakdown, and that involved the breakdown of the "engine 
configuration." (This latter modification had no effect upon the cross- 
section curves obtained. ) 


The new breakdown, as it pertained to the front portion of the tanks, 
is shown in Figure G.4-1. 


The computations were performed for: 


1, The tanks with the warhead (Configuration IV) attached, 


173 


5 [= (IRE 


S|ECIRIE 1 
UNIVERSITY OF MICHIGAN 
2428 -3-T 








Hat-Section 


R-W Re-entry Body . 


All Dimensions Shown Are in Meters 
Central Section (Cylinder: Radius =1.52, Length = 12.6) 


(See Figure G.1-1 for Breakdown of 
Aft Portion of Tank) 


FIG. G.4-1 GEOMETRY OF THE MODIFIED TANKS USED FOR COMPUTATIONAL PURPOSES 
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2. The tanks with the warhead removed, the resulting "hole’' 
being the mirror image of the rear of Configuration IV, 
and 


3. The tanks with warhead and hat section removed. 


The lengthening of the hat section led to a slightly larger value for 
the cross-sections for Cases 1 and 2 at ¥ = 80°. Otherwise, the only 
noticeable changes from the results of the Method 2 approach are in the 
cross-sections for wp = 07: 


For Case 1, the o, of Section G.3.2 was replaced by the nose con- 
tribution of Configuration IV (App. A); and for Case 2, o, was replaced 
by the wire loop formula used in Method 2 plus a contribution from the 
inside of the hole (estimated to be equivalent to that obtained from a 
sphere of radius 0.2 m, the radius of the rear hemisphere of Configura- 
tion VI). Case 3 was identical to the Case 3 of Method 2. 


The results obtained are shown in Figure G.4-2. 


175 


S 2 CIRE 


CROSS-SECTION (IN SQUARE METERS) 


S5|= CIRE 
UNIVERSITY OF MICHIGAN 
2428 -3-T 


Cross-Section of Tanks with Warhead (Conf. IV) 
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FIG. G.4-2. CROSS-SECTION OF THE TANKS AT 200 Mc (METHOD 3) 
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APPENDIX H 


IONIZATION AND RELATED TOPICS 


H. 1 INTRODUCTION 


The ionization which appears about a high-speed missile as it de- 
scends through the atmosphere could be a contributor to the radar reflec- 
tions from the missile. An analysis of this problem was begun in Ref- 
erence | and the results of that analysis are summarized below. Other 
related investigations carried out since the publication of Reference 1 
are also discussed. 


H. 2 SUMMARY OF IONIZATION STUDIES REPORTED IN REFERENCE 1 


This section contains a summary of the work on ionization produced 
by a ballistic missile performed at The University of Michigan and re- 
ported in Reference 1. This study was necessitated by the fact that, if 
missiles produced sufficient ionization, their radar cross-section might 
be significantly enhanced. Further, if they formed persistent ion trails, 
they might be confused with meteors in radar observation. 


The work was based on the fact that, to have perfect reflection from 
an ion cloud, the electron density must exceed a critical value given by 


2 


N. = a x 10° electrons/cm?, 

where fis the observing radar frequency, measured in Mc. Operation 
at 25 Mc was assumed; in practice a higher frequency would certainly be 
used, and thus a greater electron density would be required. 


The ionization was assumed to be caused by secondary collisions, in 
which atmospheric particles struck by the missile rebound elastically and 
ionization occurs in their collisions with other atmospheric particles. The 
primary collisions were treated by two opposite approximations: rigid 
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binding, where the missile atoms may be assigned infinite mass; and no 
binding, in which they are taken as free. [Ionization in the primary col- 
lision is energetically forbidden, and no other mechanisms were clearly 
established as physically significant. 


The ion distribution about the missile nose is shown to resemble a 
hyperboloid. Supracritical density is assumed within such a contour, and 
the effect of the subcritical density outside, which would reduce the radar 
cross-section, is neglected. At other points in the computation, approx- 
imations favorable to ion production are made. A missile velocity of 
10 km/sec at an altitude of 200 km was assumed. Maintenance of a con- 
tour large enough to enhance the missile's nose-on cross-section by mm?” 
was shown to require a probability > 1 for ionization in a sufficiently ener- 
getic secondary collision. It was therefore concluded that the nose-on 
cross-section of a 10 km/sec missile at that altitude would not be enhanced 
by ionization. 


It was further shown that the broadside cross-section would not be 
enhanced, since few electrons would be produced along the sides of the 
missile. Also, the electron density produced is too small to form a 
Significant trail behind the missile. 


Since the relative importance of the possible mechanisms for ioniza- 
tion is affected by altitude and speed, the fact that meteors are observed 
by radar reflection does not constitute an objection to the conclusions 
reached in Reference 1. Meteor observations by radar usually take place 
for meteor altitudes lower than 120 km (Ref. 32) and most meteor speeds 
greatly exceed the missile speeds considered. (Reference 33 contains a 
summary of the theoretical work done on the meteor problem. ) 


H.3 SHOCK WAVES 


H. 3.1 Reflections from Shock Waves 


Reflections from shock waves were also considered in Reference 1 
where it.was concluded that, for almost all missile cross-section prob- 
lems of interest, the effects of shock fronts may be neglected. 
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H. 3.2 Shape of Shock Waves During Ascent and Descent of an ICBM 


It has been conjectured that the shape of the shock wave associated 
with an ICBM might differ on the ascending and descending portions of its 
trajectory. Using trajectory data given (Ref. 34) for Atlas, this section 
will demonstrate the validity of the conjecture. 


The problem is simple for the initial portion of the ascent. Once the 
transonic region is traversed, a normal, attached, acute-angled shock 
will be present. On the basis of some photographs made at the Aberdeen 
ballistic range, and quoted by von Karman (Ref. 35), it is estimated that 
the normal pattern is established by about Mach 2. This pattern should 
persist until the missile leaves the continuum regime of aerodynamics 
and enters the slip-flow regime. The shock will then become diffuse, 
gradually reaching a state in which its existence is questionable. 


Using the Atlas trajectory data, the time history of the shock during 
the missile's ascent is determined below. 


First, consider the 5500-nautical-mile trajectory data given on page 
II-114 of Reference 34. After about 60 seconds, the transonic region has 
been traversed. A reasonable criterion for the lower limit of the slip- 
flow regime is (Ref. 17) M \/g = 0.02, where M = free-stream Mach 
number, \ = mean free path, and g is a characteristic length of the missile. 
The velocity of sound and mean free path at higher altitudes are taken 
from a Rocket Panel paper (Ref. 22). Computation then shows this missile 
to enter the slip-flow regime at an altitude of about 70 km, about 160 sec- 
onds after launching. Thus, on the way up, the missile is subsonic for 
the first 40 seconds, transonic for the next 20 or so, and finally exhibits 
Supersonic flow with a narrow shock for about the next 100 seconds. 
Starting at about that time, the shock gradually becomes diffuse and dis- 
appears. 


We consider next the ascent of the midwing-configuration glide mis- 
Sile (p. IIJ-119, Ref. 34) for 5500-nautical-mile range. Here the tran- 

sonic region is passed about 60 seconds after launching, and this config - 
uration never attains sufficient altitude to leave the continuum regime of 


aerodynamics. Thus, it retains its narrow conical shock. 
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The powered-flight trajectory for maximum range is given on page 
II-125 of Reference 34. This missile's speed exceeds Mach 2 after 60 
seconds, and it leaves the continuum regime about 140 seconds after 
launching, at an altitude of about 80 km. 


Finally, consider (p. II-127, Ref. 34) the trajectory for the 1000- 
nautical-mile range. Here, Mach 2 is passed after about 62 seconds, 
and the missile remains always in the continuum regime. Now consider 
the descent portions of these trajectories. For the glide missile, and 
the 1000-nautical-mile trajectory, unless a blunt re-entry shape gives 
rise to a detached shock, the shock configuration of the ascent is re- 
tained. This, however, is not true in the other two cases. For those 
cases, consider the possibility that a gas cap forms before the missile 
re-enters the continuum regime where a shock can exist. If this should 
happen, the shock would form around the rounded contour of the gas cap 
and therefore be rounded or possibly detached. 


The criterion for the existence of a gas cap, formed by trapped air 
molecules, is (Ref. 36) 


AV, | 
rv 


where \ is the mean free path and V the thermal velocity of air con- 
stituents, v the missile velocity, and r the missile nose radius if it 
were spherical. Since the authors have not located the assumptions 
underlying the analyses of Herlofson (Ref. 39), they have not followed 
his work here, but have used the less recent Reference 36. Here r will 
be taken as 3 inches; for an order of magnitude calculation this should 
suffice. 


For the 5500-nautical-mile powered-flight trajectory, use of the 
Rocket Panel data and the trajectories given by Convair indicate that 
the gas cap begins to form at about 110 km. The shock, however, 
doesn't begin to form until the missile has descended to about 60 km. 
Therefore, since the gas cap exists prior to the shock formation, a 
rounded or detached shock is to be expected. 
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Finally, examine the powered-flight trajectory for maximum range. 
The Convair trajectories do not present the portion of the flight imme- 
diately preceding re-entry, but the curves given are sufficiently similar 
to those for the 5500-nautical-mile range that the conclusions drawn for 
the latter retain qualitative validity. 


H:4 WORK PERFORMED AT OTHER ESTABLISHMENTS 
H.4.1 Lincoln Laboratory of the Massachusetts Institute of Technology 


Reference 1 contains a list of references containing additional ma- 
terial on ionization and related studies. Included here is the work of 
the Lincoln Laboratory of the Massachusetts Institute of Technology. 
Their work on meteors and the ICBM is reported in Reference 37. Some 
of the conclusions reached to date are: 


1, For the megawatt radars being contemplated for ICBM, 
one may expect meteor echoes at the rate of several 
thousand per hour. These meteor trail echoes will 
have a Doppler character during the trail formation, 
and will persist as relatively stationary target echoes 
for several seconds after formation. Unless the radar 
is specifically designed to reject meteor echoes, these 
can be a serious clutter and false-alarm problem. 


2. On the basis of the available meteor-rate data it ap- 
pears that the number of meteors seen by a radar 
should not increase with antenna gain but will in- 
crease with transmitter power. Hence, to get longer- 
range radars without increasing meteor background, 
one has more to gain from higher gain antennas than 
higher power transmitters. 


3, The need for further experiments is indicated, par- 
ticularly with higher-powered radars and shorter 
wavelengths than have been used until now. Such 
data are needed to establish the extent to which the 
meteor background is cut down by using higher fre- 
quencies and to determine the distribution of meteors 
corresponding to the smaller ion density trails. 
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4. The probability of a meteor hitting an offensive missile 
is negligible. 


Other work performed at the Lincoln Laboratory on the ICBM prob- 
lem is reported in Reference 38. 


H.4.2 Poulter Laboratories of the Stanford Research Institute 


Information relative to experimental work on "Extreme Velocity 
Pellets’ performed at the Poulter Laboratories of the Stanford Research 
Institute has been received in personal correspondence from Dr. Poulter 
to one of the authors. Dr. Poulter states that much work remains to be 
done to develop the controls for velocity, pellet shape, etc., but that 
they believe the solutions of these problems are rather straight-forward 
and that the methods which they have employed offer a possibility of 
obtaining any fragment velocity that will be required in missile studies. 


By using a Mach stem detonation along the axis of a comparatively 
small-angle concave conical detonation front (apparently one of the most 
effective means of collecting and concentrating energy into a very small 
space) and using 1.5 pounds of explosive, the Poulter Laboratories have 
consistently been able to photograph with the framing camera 0, 2-gram 
pellets having velocities of approximately 12 mm/p sec. In one case, 
the pellet was traveling in a lucite tube in which the air pressure was 
reduced to 50 microns. The luminous pellet, at the time the pictures 
were taken, was traveling at the tip of the jet. From the scale in the 
pictures and a camera speed of one frame per microsecond, a velocity 
of 13.2 km/sec was calculated. 


In a second case the pellet was traveling in an atmosphere of he- 
lium at a pressure of one atmosphere. The luminous conical shock 
pulse surrounding this pellet was clearly visible in the pictures and was 
seen to travel 230 mm in 17p sec, giving a velocity of 13.5 km/sec. 


One picture was taken prior to firing the charge for purposes of 


reference, the next one was taken only a few microseconds after the 
detonation front struck the rear of the metal liner and before the central 
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portion of the center of the liner had moved perceptibly. An approxi- 
mate velocity can therefore be determined for the pellet in terms of the 
velocity of the central portion of the liner and the distance that the center 
of the plate has traveled. The most conservative value that it is pos- 
sible to obtain from this picture is in excess of 60 km/sec. Thus, the 
Poulter Laboratories feel that this offers a feasible method of obtaining 


extremely high-velocity pellets in the range of from 10 to more than 
50 km/sec. 
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